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Abstract 

The recent observation by the IceCube experiment of cosmic neutrinos at energies up to a 
few PeV heralds the beginning of neutrino astronomy. At such high energies, the ‘conventional’ 
neutrino flux is suppressed and the ‘prompt’ component from charm meson decays is expected to 
become the dominant background to astrophysical neutrinos. Charm production at high energies 
is however theoretically uncertain, both since the scale uncertainties of the NLO calculation are 
large, and also because it is directly sensitive to the poorly-known gluon PDF at small-x. In this 
work we provide detailed perturbative QCD predictions for charm and bottom production in the 
forward region, and validate them by comparing with recent data from the LHCb experiment 
at 7 TeV. Finding good agreement between data and theory, we use the LHCb measurements to 
constrain the small-x gluon PDF, achieving a substantial reduction in its uncertainties. Using 
these improved PDFs, we provide predictions for charm and bottom production at LHCb at 13 
TeV, as well as for the ratio of cross-sections between 13 and 7 TeV. The same calculations are 
used to compute the energy distribution of neutrinos from charm decays in pA collisions, a key 
ingredient towards achieving a theoretically robust estimate of charm-induced backgrounds at 
neutrino telescopes. 
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1 Introduction 


The recent observation of very high-energy cosmic neutrinos by the IceCube experiment at 
the South Pole marks the beginning of neutrino astronomy [^[^. The most recent (2010-12) 
dataset contains 37 neutrino candidates with energies between 30 and 2000 TeV, and arrival 
directions consistent with isotropy. At these high energies, the ‘conventional’ atmospheric neu¬ 
trino flux, arising from the decays of pions and kaons produced by the collisions of cosmic rays 
with nuclei in the atmosphere [4]-j^ is highly suppressed due to energy loss before the decays 
occur. However charmed mesons decay almost instantaneously so at high energies, despite their 
smaller production cross-section, the so-called ‘prompt’ neutrino flux from charm decays [7-16 
becomes the dominant background to astrophysical neutrinos. The prompt flux has a harder 
spectrum than the conventional flux and is thus difficult to distinguish from the expected astro- 
physical neutrinos on this basis. 

It is therefore essential to have a reliable estimate of this prompt neutrino background. Un¬ 
fortunately, charm production at high energies is affected by substantial theoretical uncertainties 
when computed in perturbative QCD (pQCD). First of all, the small value of the charm quark 
mass (m-c), close to Aqcd, leads to a large value for as{mc), which translates into substantial 
scale uncertainties in the NLO calculation. In addition, this process probes the gluon PDF at 
very small values of x, around x ~ 10“^, where there are no direct experimental constraints and 
consequently large uncertainties 17-22 . Another source of theoretical uncertainty is the choice 
of the value of rUc itself. 

For these reasons, alternative calculations based on saturation models or non-linear evolution 
dynamics have been proposed. However, these calculations are model dependent, seldom vali¬ 
dated with collider data, and often based on outdated PDF sets. A possible alternative would 
be to use high-energy resummation for heavy quark production 
approach requires a small-x resummed PDF fit 


23 , but for consistency this 


24,25 which is currently not available. While 


there are some hints for deviations with respect to fixed-order DGLAP evolution in inclusive 
there is so far no conclusive evidence that fixed-order pQCD cannot be reli- 


HERA data 126 -28 


ably applied to the region relevant for calculations of atmospheric charm production. Therefore, 
our predictions will be based on next-to-leading order (NLO) QCD, where charm fragmentation 
is accounted for either analytically or by the matching to parton showers. 

With the above motivation, in this work we provide state-of-the-art pQCD predictions for 
charm and bottom production in the forward region. Our calculations are based both on the 
semi-analytical FONLL approach [29| , as well as the fully exclusive description of the final state 
provided by the MadGraph5_aMC@NLO 30 and POWHEG Monte Carlo programs, where the 

parton shower. As input in the calculation, we 


NLO result is matched to the PythiaS 31 


use the recent NNPDF3.0 NLO PDF set 33 


32 


modern PDF sets are used, in particular MMHT14 22 


and verify the stability of the results when other 
and CTIO 


34 . 


One central ingredient of our approach is the validation of our pQCD calculations with the 
data from the LHCb experiment on charm and bottom production in the forward region at 
7 TeV 35,36 . The LHCb measurements cover a similar kinematical range as that of charm 


production relevant to the prompt neutrino background for IceCube. For instance, an incoming 
cosmic ray with energy Fi = 100 PeV corresponds to a centre-of-mass energy of ^/s = \/2m]\fE ~ 
14 TeV. Measurements of forwardly produced heavy flavour hadrons therefore provide a perfect 
environment for testing the validity of pQCD prompt neutrino flux predictions. As we will 
show, both the analytical FONLL calculation and the exclusive MadGraph5_aMC@NLO and 
POWHEG results are consistent with the LHCb charm and bottom data within theoretical 
uncertainties. Therefore, we can be confident that these calculations can be reliably applied to 
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predictions of the atmospheric prompt neutrino flux. 

The compatibility between the NLO QCD predictions and the 7 TeV charm production 
data from LHCb indicates that it is possible to use this process to constrain the small-x gluon 
PDF 37 . To partially cancel the large scale uncertainties of the NLO calculation, we con¬ 


struct normalised differential cross-sections using a fixed bin as reference. We then include the 
LHCb charm data into NNPDF3.0 fit using the Bayesian reweighting method [M 40 , find¬ 
ing a substantial reduction of the small-x gluon PDF uncertainties. The resulting PDF set, 
NNPDF3.0-|-LHCb, is particularly suitable for providing predictions for both heavy quark pro¬ 
duction within the LHCb acceptance at 13 TeV, as well as to provide a reliable estimate of the 
rate of high energy neutrino production in pA collisions relevant for estimations of the prompt 
neutrino flux at IceCube. 

In this work we also provide detailed predictions for charm and bottom production at LHCb 
Run H, using the improved NNPDF3.0-|-LHCb PDFs, including the evaluation of theoretical 
uncertainties arising from missing higher-orders, PDFs, and the value of the heavy quark mass. 
Our results are tabulated using the binning scheme adopted for the 7 TeV measurements, and 
predictions for other binning choices are available upon request. In addition, we also provide 
predictions for the ratio of differential distributions of charm and bottom production between 
13 and 7 TeV, R 13 / 7 , which provides complementary information on PDF discrimination 


41 


After computing this observable and its corresponding theoretical uncertainty for B and D 
mesons, we apply our calculations to the LHCb 7 TeV data to provide robust predictions for 
the fiducial cross-sections within the LHCb acceptance for Run H. These predictions are useful 
for estimating B and D yields at 13 TeV, which can in turn be used to assess the statistical 
precision of future measurements — such as rare B decays for example. 

Using the same theoretical set-up as outlined for the LHC calculations, we provide predictions 
for the neutrino energy spectrum arising from the decays of charmed mesons in high-energy 
proton-air collisions. These results are an important ingredient for the computation of the 
expected number of prompt neutrino events at IceCube. While it is beyond the scope of this 
paper to compare with the IceCube measurements, our pA —)• vX cross-sections are available 
in the form of an interpolation code for the relevant range of incoming cosmic ray energies. 
These results can be used as an input for well-established frameworks such as the Z-moment 


approach 10,13 to construct predictions for IceCube. 

The outline of this paper is as follows. In Sect. we review the framework for pQCD 
computations of heavy quark production and perform an extensive comparison with the LHCb 
7 TeV data on charm and bottom cross-sections. In Sect. [3] we include the normalised LHCb 
charm data into NNPDF3.0 using the Bayesian reweighting method, obtaining an improved 
PDF set with reduced uncertainties in the small-x region which will be the central ingredient 
of our subsequent calculations. In Sect. we provide predictions for heavy quark production 
within the LHCb acceptance at 13 TeV, as well as ratio of 13 over 7 TeV cross-sections. In 
Sect. we present our predictions for the energy distributions of neutrinos from charm decays 
in pA collisions for a range of incoming cosmic ray energies relevant for neutrino telescopes. 
In Sect. we summarise our findings and discuss possible next steps. Appendix contains a 
tabulation our theory predictions for charm and bottom production at LHCb at 13 TeV, as well 
as the ratio of cross-sections between 13 and 7 TeV. 


2 Heavy quark production in the forward region and LHCb data 

In pQCD, the NLO calculation of heavy quark pair production in hadronic collisions has been 
available for a long time, both at the level of total inclusive cross-sections |42|, and of differential 
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distributions (43]-l4^. Subsequently, the fixed-order calculation has been improved with the 


resummation of soft gluons at NLL 47,48 and NNLL |49l 50] accuracy. Another way of refining 


the fixed-order result is by matching it to the massless calculation, valid in the limit where 
the heavy quark transverse momentum (p^) greatly exceeds the heavy quark mass (m/j), thus 

and has the 


29 


51 


52 


obtaining a result which is valid both at small and at large values of 
benefit of reduced of scale uncertainties as compared to the NLO calculation. More recently, 
the next-to-NLO (NNLO) calculation for inclusive heavy quark pair production has become 


available 53 -^, and results for the differential distributions for the case of top quark production 
have also been presented 56,57]. These calculations will eventually be applied to charm and 


bottom production as well. 

In this section, we begin by discussing our set-up for providing pQCD calculations of charm 
and beauty production, and their subsequent fragmentation and decay. We then demonstrate 
that the kinematic coverage of charm production at LHCb data overlaps with that relevant for 
the calculation of prompt neutrino fluxes at IceCube. With this in mind, we present a detailed 
comparison of the pQCD calculations for charm and bottom production in the forward region 
with the 7 TeV LHCb data, and examine relevant sources of theoretical uncertainty. Throughout 
this work, the NNPDF3.0 NLO PDF set will be used as a baseline for our predictions, and we 
also study the dependence of our predictions on the choice of input PDF set. 


2.1 Heavy quark production in the forward region 

In this work we will provide pQCD predictions of heavy quark pair production using three 
different approaches: FONLL, POWHEG and MADGRAPH5_aMC@NLO. We discuss briefly each 
of these approaches in turn. A similar comparison between different calculations for heavy quark 
production at the LHC and other hadron colliders, focussed on data in the central rapidity region, 
was presented in . 


FONLL 29 51,59 


fixed-order calculation 43 


is a semi-analytical calculation based on the matching of the NLO 
including full dependence on the heavy quark mass mh, with 


the resummed NLL calculation where the heavy quark is treated as a massless parton. 
This matching allows a consistent description of the spectrum, from low to high trans¬ 
verse momenta0 The fragmentation of heavy quarks into heavy flavored hadrons is then 
described analytically 61 , with parameters extracted from LEP data. It is also possible 


to include the decays of the D mesons using this approach. 

In the region relevant for the LHCb data, where does not greatly exceed mh, the FONLL 
result corresponds to the fixed-order NLO massive calculation, and thus for simplicity in 
this work by “FONLL calculation” we denote the fixed-order NLO obtained from the 
FONLL code. 


The POWHEG Ri2-64 


method allows NLO calculations to be matched to a Monte Carlo 
parton shower. In the case of heavy quark production |65|, the massive NLO calculation 


performed in a fixed-flavour scheme is matched achieving NLO-I-LL accuracy — thanks to 
the resummation achieved by the parton shower. The fragmentation and hadronisation of 
heavy quarks into heavy hadrons and their subsequent decay into leptons is then modeled 
by the specific parton shower which has been matched too, with modelling parameters 
tuned to data. In this work we use POWHEG matched to the PythiaS shower 
using the Monash 2013 tune for the modelling of the soft and semi-hard physics 
will refer to this set-up as the POWHEG calculation. 


31,32 


We 


The FONLL approach can also be applied to other processes, such as DIS structure functions 


60 


5 



























MADGRAPH5_aMC@NLO [30] provides automated calculations of arbitrary processes at 
LO and NLO, both at fixed-order and matched to a variety of parton showers using the 
MC@NLO method 67 . For consistency with the POWHEG calcnlation, the PythiaS 


parton shower and Monash 2013 tnne are also used for this prediction, therefore treating 
charm hadronisation and decay with universal settings. Note that the MC@NLO and 
POWHEG methods to match fixed-order calculations with parton showers are different, 
and thus the spread between the two calculations provide an estimate of the underlying 
theoretical nncertainties introdnced by the various matching processes. This set-up is 
referred to as the aMC@NLO calculation. 


In the kinematic region relevant for charm and bottom prodnction at LHCb, the effects 
of parton shower resummation in POWHEG and aMG@NLO are expected to be moder¬ 
ate, and thus the comparison of the three generators allows a meaningful validation of 
the pQCD calculations for the heavy qnark production and fragmentation nsing three 
independent approaches. 


The following common set of theory input parameters are adopted for all three calculations: 


As the input set of parton distribntions, we use the NNPDF3.0 NLO PDF set with 
five active flavours (n/ = 5). The dependence of the results with respect to the choice of 
input PDF set will be discussed in Sect. 2.4 and comparisons with recent PDF fits will 
be made in Sect. At the LHC, charm and bottom pairs are predominantly produced 
through the gluon-gluon initial state, and therefore our calculation will be sensitive to the 
details of the gluon PDF and the associated uncertainties at small-x. 


Charm prodnction in the presented FONLL predictions only includes the matching between 
the Uf = 3 to the n/ = 4 schemes, so in principle one shonld use a n/ = 4 PDF set for 
consistency. However, it has been verified that the results are unchanged in the latter 
scenario: differences between using FONLL with n/ = 4 and n/ = 5 PDFs for charm 
production are at most 1.5% at the highest values of p!^ covered by the LHCb data, much 
smaller than any other theoretical or experimental uncertainty. 


In the case of both POWHEG and aMC@NLO calculations, the matching between schemes 
is not included. We have verified however, by explicitly including these terms in the 
calculation, that such effects are also in this case unimportant. In partic- 


POWHEG 71 


nlar, the effect of including the nj = 3tonj = 5 compensation terms in the POWHEG 
calcaulation with a n/ = 5 PDF set leads to an increase in scale variation of (2-3)% above 
mb, while the central valne is essentially unaltered (< 1%) due to a compensation of Uf 
dependent as modifications and a depletion of the gluon PDF due to g —)• QQ splittings. 


For completeness, we provide here the explicit expressions of the compensation terms 
that must be dynamically applied to the partonic heavy-quark production cross-section to 
transform from the nj to the Uf + 1 scheme: 


( 0 ) ‘^'^FOisip'ji) 


Log 



- ( 0 ) ‘^'^POisip'ji) 
""39 3vr 




( 1 ) 


These expressions are valid for the choice pp and pR > nrg. If only the value of pp exceeds 
mg, then only the /ii?-dependent correction to the gluon-gluon induced process should be 
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applied (and similarly for the /i/j-dependent corrections). In the case of charm production, 
li pLp and ^lr exceed then the corrections 0 should be applied at both charm and 
bottom thresholds. 

In addition, let us recall that differences between nj = 4 and n/ = 5 PDFs are only 
sizeable far above the bottom threshold [68| , thus not relevant for the analysis of the 
LHCb production data. 

• The value of the strong coupling constant is taken to be as{mz) = 0.118, consistent with 
the latest PDG average [^. The uncertainties due to the uncertainty of the value of 
as{mz) are negligible as compared to other sources of theory uncertainty and are thus not 
considered here. 

• Concerning the treatment of as{Q), in this work we always use consistently the same heavy 

flavour scheme as the corresponding input PDF set. Since we use nj = 5 PDF sets, then 
as{Q) runs with up to n/ = 5 active flavours depending on the value of Q. Close to the 
charm threshold, as ■' {Q) and a^ ■' (Q) are extremely similar by construction. 

Note also that the VFN running of as{Q) is essential to obtain agreement with the PDC 
global average of as{mz)- using the nj = 3 scheme all the way up to Q = mz will lead to 
a value of as{mz) much smaller than the PDC average. 

• The central renormalisation and factorisation scales are varied event-by-event, and taken 
to be 



To estimate the size of missing higher-order corrections, and ^r are varied by a factor 
of two around the central scale, with the restriction \l‘l< [ipjl^R < 2 to avoid introducing 
artificially large logarithms. Uncertainties computed in this way are referred to as scale 
uncertainties. 


The charm quark pole mass is taken to be rric = 1.5 ±0.2 CeV, while for the bottom quark 
pole mass we use mb = 4.75 ± 0.25 CeV. The uncertainty of me and mb will be included 
in the theory uncertainty of our calculation. While it should be possible to reduce the 
theory uncertainty due to the choice of heavy quark masses by using calculations in the 
MS scheme 70 , where the latest PDC values are mdme) = 1.275 ± 0.025 CeV and 
iRbiiT^b) = 4.18 ± 0.03 CeV [^, this would not affect our results since the uncertainties 
due to 6mc (and even more due to 5mb) are subleading as compared to other theory 
uncertainties. 


The fragmentation probabilities /(c 

taken to be the same as those of the LHCb measurement 
/(c ^ D±) = 0.246, /(c ^ Df ) = 0.080, and f{c Ac) = 0.094. 


36 


VIZ. 


D) for the different types of charmed mesons are 

f{c D°) = 0.565, 
When uncertainties 

are considered, the sum of these fragmentation probabilities is consistent with unity. In 
comparison to the other sources of theoretical uncertainty, the impact of the uncertainty 
of these values for the considered observables is negligible. 


• When semi-leptonic decays of D hadrons are considered, the following branching fractions 
are enforced: B{D^ —)■ viX) = 0.101, B{D^ —)■ RiX) = 0.153, B{Df —)• uiX) = 0.06, and 
H(Ac —7- viX) = 0.02. Combined with the fragmentation probabilities, this corresponds to 
a partial decay width F(c —)• i>iX)/T{c —)• anything) = 0.102 for prompt D hadron decays. 
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Figure 1: Contour plot for the values of {xi,X 2 ) sampled in the LO calculation of charm (upper plots) 
and bottom (lower plots) production at 7 TeV, within the LHCb fiducial acceptance. The calculation 
has been performed with POWHEG using the NNPDF3.0 LO set. The regions in red indicate where 
the PDFs are sampled more frequently, while those in blue indicate less frequent sampling. The left plots 
have been computed in the full fiducial region, while the right plots are restricted to the forward region 
4.0 < y < 4.5. 


• The fragmentation probabilities f{b —)■ B) for bottom mesons are taken to be f{b —)• Bu) = 
f{b —)■ Bd) = 0.337, as determined by the LHCb analysis of Ref. (7^ . 

2.2 Sensitivity to the small-x gluon PDF 

In order to better understand the relation between heavy quark production kinematics and the 
gluon PDF, it is useful to determine the coverage in the {xi,X 2 ) plane of the LHCb charm and 
bottom measurements, where xi and X 2 are the values of Bjorken-a: corresponding to the PDFs in 
each of the two incoming protons. This coverage is illustrated by the various contour plots shown 
in Fig. These plots contain the values of {xi,X 2 ) sampled by the LO calculation of charm 
(upper) and bottom (lower) production at 7 TeV, within the LHCb acceptance. In the left plots, 
and B^ hadrons are required to be within the LHCb rapidity acceptance (2.0 < y < 4.5) 
and have been restricted to a low px region [px < 8 GeV). In the right plots, the hadrons are 
further restricted in rapidity to the most forward region with 4.0 < y < 4.5. The calculation 
has been performed with POWHEG using NNPDF3.0 LO. In all plots, the contours have been 
normalised to the corresponding fiducial region, and therefore the regions in red indicate where 
the PDFs are sampled more frequently, while those in blue indicate less frequent sampling. Note 



















































































that due to the asymmetric acceptance of LHCb, events with xi > X 2 , where the first parton is 
a constituent of the proton travelling in the direction of the LHCb detector (positive rapidity), 
will be typically selected. 

As shown in Fig. measurements of charm production probe average values of Bjorken-x 
as low as (X 2 ) — 4.6 • 10“®, and even knowledge of the gluon PDF for values below x < 10“® is 
required for particular bins. This is demonstrated by the plot restricted to the forward region, 
where (X 2 ) — 1.5 • 10“^. In this region, there is very limited direct experimental information, 
since HERA inclusive structure function data |26| is only available down to x 

min 6-10 For 

this reason, it is of paramount importance to validate our pQCD calculation with the LHCb 
data itself, since we are using as input PDFs in a region where uncertainties are extremely large. 
In contrast, the situation for bottom production is under better control since (X 2 ) — 1-3 • 10“'^, a 
region well covered by the HERA data. This said, for bottom production in the most forward bin, 
4.0 < y < 4.5, we find that (X 2 ) — 4.7 • 10“®, just below the limit of HERA data, demonstrating 
that PDE uncertainties also have a sizable impact in this region. 

To better illustrate this point, and bearing in mind that heavy quark production at the LHC 
is driven by the gg luminosity, it is useful to quantify the PDE uncertainties of the NNPDE3.0 
gluon, and compare this to other NLO PDE sets. To ease these comparisons, we use the APFEL 
Web on-line PDE plotter [^[^. In Eig. i we show a comparison of the gluon PDEs evolved 
to the scale Q = 1.4 GeV (corresponding to a typical value of the charm mass) between the 
NNPDE3.0 and (from top to bottom) the CTIO 


20 and MMHT14 22 


NLO PDE sets. In 

each case, the bands correspond to the 68% confidence level for the PDE uncertainties. The 
right plots of Eig. show the same comparisons now performed at the scale Q = 4.5 GeV, a 
value typical of the bottom quark mass, shown as ratios with respect to the central NNPDE3.0 
prediction. 

As can be seen, in the region relevant for charm production at LHCb, with X 2 < (X 2 ) — 
4.6 • 10“®, the gluon PDE uncertainties are extremely large. On the other hand, for the region 
relevant for bottom production, with X 2 < (X 2 ) — 1.3 • 10“^, PDE uncertainties are moderate, 
thanks to the constraints from HERA data. Importantly, as shown in Eig. the description 
of the gluon PDE at small-x is quite similar, both in terms of the central value and associated 
uncertainty — particularly for the comparison between NNPDE3.0 and MMHT sets. As will be 
shown explicitly, this agreement implies that predictions for charm and bottom production at 
LHCb obtained with NNPDE3.0 will be similar to those obtained with CTIO or MMHT14 as 
input PDE sets. 

In Eig. 1^ we show the comparison between the n/ = 4 and n/ = 5 gluon and up quark 
NNPDE3.0 NLO PDEs as a function of Q, for a reference value x = 2 • 10“®, in the kinematical 
region relevant for charm production at LHCb. We see that the differences between the nj = 4 
and Uf = 5 schemes are much smaller than the associated PDE uncertainties. We have also 
explicitly verified that either using n/ = 4 PDEs in the FONLL calculations or including the 
Uf —7- Uf + 2 scheme transformation terms in POWHFG leads to negligible modifications of our 
results. These considerations justify our choice of the NNPDF3.0 NLO Uf = 5 set as baseline 
in our calculations. 

The fact that gluon PDF uncertainties in the region relevant for charm production at LHCb 
are large indicates that these measurements can be used to provide information on the poorly 
known small-x gluon. This constraining potential has been recently verified by the PROSA 

In Sect. we will study the impact of 


analysis |37| based on the HERAfitter framework 75 


the LHCb charm data in the NNPDF3.0 NLO global analysis using the Bayesian reweighting 
method. 
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Figure 2; Left plots: comparison of the small-x gluon PDFs at Q = 1.4 GeV between NNPDF3.0 and 
(from top to bottom) CTIO and MMHT14. PDFs are compared in an absolute scale, and the bands 
indicate the PDF uncertainties. Right plots: the same comparisons performed at Q = 4.5 GeV, now 
shown as ratios with respect to the central NNPDF3.0 prediction. 


2.3 Comparison with the LHCb data 

We now perform a detailed comparison of the pQCD calculations of charm and bottom produc¬ 
tion in the forward region with the most recent LHCb data 
performed at the level of double differential distributions. 


35 36 . The comparisons will be 


d?a{D){y,PT) 

dy^dpS 


and 


d^a{B){y,pT) 

dy^dpS 


( 3 ) 


For all mesons, we have also checked that good agreement is obtained for the total cross-sections 
in the fiducial region. 

For D mesons, we restrict the comparison to the case of the higher-statistics final states, 
namely and while for the beauty mesons we will show results only for production. 
For each calculation, we provide the central prediction as well as the contribution arising from 


the various sources of theoretical uncertainty as outlined in Sect. 2.1 


The comparison between the FONLL calculation and the LHCb charm production data is 
shown in Fig. We show the results for the most central bin, 2.0 < y < 2.5 and a forward bin, 
3.5 < y < 4.0, both for the and the measurements. In Fig. statistical and systematic 
uncertainties have been added in quadrature for the experimental data, while for the theory 
uncertainties we show both the scale uncertainty alone and also the sum in quadrature of scale 
and PDF uncertainties. 
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Figure 3: Comparison between the nj = 4 and Uf = 5 PDFs from the NNPDF3.0 NLO set, as a function 
of Q for a; = 5 • 10“^, in the region relevant for forward charm production at LHCb. We show the gluon 
(left plot) and the up quark (right plot), normalized to the central value of the Uf = 5 set. 


The agreement, within uncertainties, between the LHCb data and the NLO pQCD prediction 
across the entire kinematic range demonstrates the applicability of this approach to forward 
charm production. The total theoretical uncertainty is dominated by scale variation, except in 
the low pt where the large gluon PDF uncertainty at small-x becomes comparable to the scale 
variation or even dominant. Similar satisfactory agreement is found for the other data bins not 
shown in Fig. 

Given the compatibility of the charm production data and theory prediction provided by 
FONLL, we now compare these predictions to those obtained with the NLO Monte Carlo ap¬ 
proaches, aMCONLO and POWHEG. First of all we compare the FONLL results with the 
aMC@NLO calculation. For simplicity, we only provide results for mesons. The comparison 
is shown in Fig. clearly, there is good agreement between the central values of the two cal¬ 
culations. For the total theory uncertainty band there is also reasonable agreement, with the 
aMC@NLO band being typically larger than, but still consistent, with the FONLL result. In this 
comparison the theory uncertainty band is obtained from adding scale and PDF uncertainties 
in quadrature. The corresponding comparison between the two MC generators, aMC@NLO and 
POWHEG, is shown in the lower plots of Fig. Reasonable agreement is also found between 
the central predictions, well within the uncertainty bands. We note that scale uncertainties tend 
to be slightly larger in the POWHEG calculation^ 

In Fig. we perform the same comparison between the three calculations as shown in Fig. 
but now normalising each prediction to the corresponding central value. This way we can 
gauge how the total theory uncertainty band compares among the three calculations. The 
total uncertainty is similar for POWHEG and aMG@NLO calculations. Notably, the scale 
uncertainties of the POWHEG and aMG@NLO calculations tend to be larger than those of 
FONLL, especially in the upper variations in the moderate and high pT region. While the origin 
of these differences remains to be understood, it might be related to the fact that FONLL is 
a fixed-order calculation while POWHEG and aMG@NLO are matched to parton showers, and 
this matching may induce additional theoretical uncertainties. Indeed, we have verified that the 


^This has been traced back to a different solution of the RG equations for the running of as{Q) used in the 
POWHEG calculation, leading to formally subleading corrections which are numerically important at Q ~ rric- 
As opposed to aMC@NLO and FONLL, where as{Q) is consistently extracted from the PDF set that is being 
used via the LHAPDF6 76 interface, POWHEG uses its own internal routine for the running of q:s(Q). We 
thank Emanuele Re for clarifications about this point. 
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Figure 4: Comparison between the LHCb data on D meson production and the FONLL calculation 
using NNPDF3.0 as input. We show the results for the most central bin, 2.0 < y < 2.5 (left column) and 
a forward bin, 3.5 < y < 4.0 (right column), both for data (upper row) and the data (lower row). 
The solid error band is obtained from the sum in quadrature of PDF and scale uncertainties, while the 
hatched band is only the scale variation component. 


scale uncertainties of the fixed-order NLO computation of differential cc production (without 
fragmentation) in aMC@NLO reproduces those of FONLL to a few percent. 

From Fig. we see that the FONLL semi-analytical calculation exhibits smaller theoretical 
uncertainty, and for this reason, in the following Section we will use the FONLL predictions to 
quantify the constraints of the LHCb charm production data on the NNPDF3.0 small-x gluon 
PDF. 


We now begin the comparison between the LHCb data and the various theoretical calcula¬ 
tions for the case of B meson production. For simplicity, we show results only for mesons, 
though similar agreement has been found for the other B mesons. As compared to the case of 
the D mesons, we expect a reduction of the theory uncertainties for several reasons: the calcu¬ 
lation is performed at a higher scale + Pj^b^ compared to the charm production case, 


leading to an improved convergence of the perturbative expansion; the relative 
uncertainty of the value of mb is smaller; and larger values of xi ^2 are probed within the proton, 
a region well covered by HERA data as illustrated in Fig. [^and Fig. 

In Fig.j^we show the comparison of the LHCb data for B^ meson production, both for central 
and for forward rapidities, with the corresponding POWHEG and aMC@NLO calculations. The 
indicated theory uncertainty band includes only the scale uncertainties, and we have verified 
that PDF uncertainties are not so relevant in this case. As in the case of charm, satisfactory 
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Figure 5: Comparison between the FONLL and aMC@NLO (upper plots) and between the POWHEG 
and aMC@NLO (lower plots) calculations for production in the same kinematics as the LHCb data 
of Fig. 1^ using NNPDF3.0 NLO. The total theory uncertainly band is obtained by the addition in 
quadrature of scale and PDF uncertainties. 


agreement between theory and data for B meson production in the forward region is found. 
There is also a substantial reduction of the theory uncertainty as compared to the D meson 
case. The POWHEG and aMG@NLO predictions are in reasonable agreement within the theory 
uncertainty band. 

To better assess the differences between the two NLO matched calculations, we compare 
them again in Fig. this time with the distributions normalised to the central POWHEG 
prediction. The aMG@NLO and POWHEG predictions agree across the considered kinematic 
range, with the POWHEG prediction favouring a slightly larger cross section in the low pT 
range. In comparison to the charm results, Fig.[^ the reduction of scale uncertainties is evident, 
since now the scale variation amounts to an uncertainty of ~ 40%. We can conclude that the 
pQCD description of B meson production in the forward region is completely satisfactory, and 
that theory uncertainties are substantially reduced as compared to charm production. 

In this section we have restricted our study to 7 TeV, the only centre-of-mass energy for which 
LHCb measurements are currently available. Predictions for double differential distributions at 
13 TeV, as well as for the ratio of cross-sections at computed at 13 over 7 TeV, will be provided 
in Sect. 14.11 and 14.21 

2.4 PDF dependence of heavy quark production at LHCb 

The results shown so far in this Section have been computed using the NNPDF3.0 NLO set. We 
have verihed that the pQCD predictions for heavy quark production are affected by a sizeable 
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Figure 6: Comparison between the total theoretical uncertainty (sum in quadrature of scale and PDF 
uncertainties) for the kinematics of production at LHCb. The results for the three calculations, 
aMC@NLO, POWHEG, and FONLL calculations, are normalised to the respective central values. 

B“ mesons, 2.0 < y < 2.5 B“ mesons, 3.5 < y < 4.0 




Figure 7: Comparison of the LHCb data on meson production, both for central and for forward 
rapidities, with the theoretical predictions from POWHEG and aMG@NLO. The theory uncertainty 
includes only scale uncertainties. 


PDF uncertainty, which arises in turn from poor knowledge of the small-x gluon PDF due to a 
lack of direct experimental constraints. In this section we study the dependence of our predictions 
on the choice of input PDF set, in particular we compare those of the baseline NNPDF3.0 to 
CTIO and MMHT14 NLO sets. The comparison of the small-x gluon PDF between these three 
sets shown in Fig. [^indicates that predictions for charm production cross-sections are expected 
to be reasonably similar. 

In Fig. [^we show the comparison of the theoretical predictions for charm production at 7 
TeV within the LHCb acceptance found using the POWHEG calculation with NNPDF3.0, CTIO 
and MMHT14 PDFs. The uncertainty band corresponds to the 68% confidence level for each 
PDF set, and the shown results have been normalised to the central value of the NNPDF3.0 
prediction. From this comparison, we see that the dependence of the charm cross-section on the 
choice of input PDF set is moderate, with the three central values consistent within large PDF 
uncertainties. Recall that at fixed rapidity, smaller values of the D meson pT correspond to 
probing smaller x values for the gluon PDF, and that, likewise, for a fixed value of pT, forward 
rapidities corresponds to smaller x values. It is therefore reasonable that PDF uncertainties are 
largest at small px and forward rapidities, as shown in Fig. 
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Figure 8: Comparison of the theoretical predictions for meson production at the LHCb kinematics 
between POWHEG and aMC@NLO, shown in Fig.[^ but now where calculations are normalised to the 
central value of the POWHEG prediction. 
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Figure 9: Comparison of the theoretical predictions for meson production at y/s=7 TeV in the LHCb 
kinematics with POWHEG, for three different NLO sets of PDFs, NNPDF3.0, CTIO and MMHT14. 
The band corresponds to the respective one-sigma PDF uncertainty for each set. Results are shown 
normalised to the central value of the NNPDF3.0 prediction. 


Even though predictions suffer from large PDF uncertainties, the central value of these 
three PDF sets are reasonably consistent. This agreement can in part be explained by the fact 
that at small-x PDF constraints in the three sets come from the same dataset, the combined 
HERA-I measurements |26] . We note that the relative size of the PDE uncertainties is similar 
for NNPDF3.0 and CTIO, while the MMHT14 uncertainty is about a factor of two smaller. 
Another feature of these predictions is the preference for the CTIO and MMHT14 central values 
towards relatively smaller and larger differential cross sections for small pT values, respectively. 
This can be traced to the relatively softer and harder gluon PDF at small-x preferred by the 
CTIO and MMHT14 respectively as compared to NNPDF3.0 — see Fig. 

We conclude from Fig. [^that although there is some dependence on the choice of input PDF 
set, these differences are small within the large intrinsic PDF uncertainties, and therefore it is 
sufficient to use a single PDF set, NNPDF3.0, as baseline in our calculations. 
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3 Constraints on the small-x gluon PDF from forward charm 
production data 

As demonstrated in the previous section, the production of charmed hadrons in the forward 
region and the associated theoretical uncertainty depends on the description of the gluon PDF 
at small-x. We now use the charm production data from LHCb to substantially reduce the small- 
X gluon PDF uncertainties. This will allow a more reliable prediction for both forward charm 
production at the LHC Run II and the prompt neutrino cross section arising from high energy 
cosmic rays — an important input for calculating the background neutrino flux at IceCube. 

The basic idea is similar to the study performed by the PROSA Collaboration , where the 
impact of forward B and D LHCb data on the low-x PDFs is studiecQ The PROSA study is 
based on the HERAfitter framework [^, and quantifies the error reduction in a HERA-only 
PDF fit when the LHCb B and D meson production data is included using the MNR code |46j in 
a FFN Nf = 3 scheme. Similarly as will be done here, theoretical uncertainties can be reduced 
by suitable normalisations. This said, there are important methodological differences in the two 
analysis (global fit versus HERA-only fit, theory calculations, data normalisation strategies), 
and so the two approaches complement one another. 

The starting point is the NNPDF3.0 NLO set, with as{mz) = 0.118, supplemented by the 
LHCb measurements of the 7 TeV differential distributions for and production 36 . The 


LHCb data will be added to the NNPDF3.0 global dataset by means of the Bayesian reweighting 
technique 39,40 . This method allows to quantify the impact of new data in a set of Monte 


Carlo PDFs without the need of redoing the full global QCD analysis, and has been used before 
in a number of related applications in order to quantify the impact on PDF fits from data for 


isolated photon production 77,78 , top quark pair production 79 , and polarised and jet 


production . As an alternative to the reweighting, it should also have been possible to use the 
aMCfast |81| interface to construct an APPLgrid 82 fast implementation of the aMC@NLO 
calculations presented in the previous section. 

As input to the reweighting, we consider the (y,pr) double differential distributions for 
and production at LHCb, but exclude the data from other final states such as D*^ and 
which are affected by larger experimental uncertainties, and therefore have reduced impact on 
the fit. These data cover a range in rapidity of [2.0,4.5] and in px of [0, 8] GeV. In total, we are 
adding A^dat = 75 new data points into the NNPDF3.0 analysis. 

For the theoretical calculations, we use the FONLL predictions, with the settings discussed in 
the previous Section. In Fig. [T^we compare the LHCb charm production data and the FONLL 
prediction for the and data. Results are shown normalised to the central value of the 
respective experimental data point. The experimental statistical and systematic uncertainties 
have been added in quadrature, and both scale and PDF uncertainties are independently shown 
for the FONLL theoretical prediction. 

It is clear by inspection of Fig. that the scale uncertainties in the NLO calculation are 
large, by as much as a factor of two for some bins. In general, they are reduced when going 
towards higher pT bins, thanks to the improved convergence of the perturbative expansion 
in this region. Although PDF uncertainties are also large, especially at low pT and forward 
rapidities where the small-x gluon is being probed, they are sub-dominant as compared to the 
scale uncertainties. This is concerning from the point of view of a PDF analysis, in which a 
scale choice for the central value of the theory prediction must be made. 


^We would like to stress that preliminary results for our work were presented already in February 2015, 
http://benasque.org/20151hc/talks_contr/179_BenasqueGauld.pdf, before the publication of the PROSA pa¬ 
per. Preliminary results of the PROSA study were also presented in 381. 


16 












7 TeV D° unnormalized 


7 TeV unnormalized 




Figure 10: Comparison between the LHCb charm production data and the FONLL calculation with 
NNPDF3.0 NLO at the level of unnormalized (absolute) cross-sections. The left plot shows the 11° data 
while the right plot corresponds to the data. Results are shown normalised to the central value of 
the LHCb data. For the FONLL calculation we show separately the scale and the PDF uncertainties. 
The data is ordered in increasing rapidity bins, an within each of these, in increasing pr bins. 


7 TeV d“ normalized 7 TeV □*' normalized 




Figure 11: Same as Fig. 10 where now both data and theory have been respectively normalised with 
respect to the bin with 7 GeV < pt < 8 GeV and 2.0 < y < 2.5 (the bin with data point index 8). 


To bypass this problem, the strategy that will be adopted in this work is to normalise all the 
data bins to that with highest p!^, [7,8] GeV, and central rapidity , [2.0, 2.5]. The rationale 
for this choice is that scale uncertainties will partially cancel in the ratio, while the cancellation 
of PDF uncertainties will not be as severe, given that different bins in probe different 

values of (x, Q^) of the gluon PDF. The reference bin has been chosen precisely for this reason, as 
PDF uncertainties for this particular bin are the smallest. Note that this is strategy is different 
as compared to the PROSA analysis 37 , where, separately for each bin in the rapidity bin 


3.0 < < 3.5 was used to normalize the data and the theory calculations. 


In Fig. [m we provide the same comparison of Fig. but this time at the level of nor¬ 
malised distributions. In Fig. im we have added in quadrature the experimental uncertainties 
in the numerator and the denominator, this being the only option since the full experimental 
covariance matrix with the information of correlations between bins is not available. Theoretical 
uncertainties are taken to be fully correlated among all the data bins. 

The comparison between Fig. 10 and Fig. |11| illustrates how after the normalisation proce¬ 
dure has been applied, scale uncertainties are substantially reduced in the low-pT and large-y 
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bins. Importantly, the PDF uncertainties are now larger than the corresponding scale and ex¬ 
perimental uncertainties in these bins, which justifies the inclusion of the normalised charm 
production cross-sections into NNPDF3.0 using Bayesian reweighting. In this respect, after the 
normalisation, the theoretical status of forward charm production becomes similar to that of 
other hadronic processes routinely included in global NLO fits, such as jet production. 

The results of the reweighting are summarised in Table The breakdown of the per 
data point of the and data before and after reweighting, as well as the number of 
effective replicas left out of the original iVrep = 100 replicas, is provided. The description of 
the normalised LHCb charm data turns out to be excellent even using the original NNPDF3.0 
set, with a value of x^/-^dat = 1.10. This is certainly reassuring, since it shows that both 
NNPDF3.0 and the FONLL calculation provide a good description of charm production in the 
LHCb acceptance. Once the data is included by the reweighting, the x?w/-^dat = 0.74 is even 
better, and the effective number of replicas is iVeg = 50, confirming that this data is indeed 
very constraining on the small-x gluon PDF. Note that since we are neglecting the correlations 
between systematics, we are underestimating the impact of these data. Future measurements 
with the full systematic breakdown should be even more powerful. 



NNPDE3.0 

NNPDE3.0-kLHCb data 

X^/A^dat for + c.c. 

1.13 

1.05 

X^/Adat for D+ -k c.c. 

1.06 

0.40 

xVAdat for 

1.10 

0.74 

Neff for 

- 

50 


Table 1: Results of the reweighting of NNPDF3.0 with the LHCb charm production data. We give the 
value of the x^/-^dat, both for the original NNPDF3.0 set, and for the reweighted NNPDF3.0-|-LHCb set, 
as well as the effective number of replicas left, iVeg. 


The impact of the LHCb charm production data into the small-x gluon PDF can be seen 
in Fig. [T^ We show the NNPDF3.0 small-x gluon, evaluated at Q = 2 GeV, compared with 
the new gluon obtained after the inclusion in the fit of the normalised LHCb charm data. As a 
cross-check, we have also verified that it is possible to unweight the results to produce a stand¬ 
alone LHAPDF6 grid for the combined NNPDF3.0-|-LHCb fit (indicated as “(unw)” in the 
plot legend). In Fig. 12 we also compare the percentage PDF uncertainties for the NNPDF3.0 
gluon with and without the inclusion of the LHCb data, which quantify the reduction of PDF 
uncertainties at small-x. 

We see that the impact of LHCb data is negligible for x > 10““^, where most of the HERA 
data is available, but becomes substantial for x < 10“^, where the previously large PDF un¬ 
certainties are dramatically reduced. For instance, for x ~ 10“®, the PDF uncertainties in the 
gluon PDF are reduced by more than a factor three. We also note that the central value at 
small-x of the gluon PDF preferred by the LHCb charm data is less steep than that of the global 
fit, although fully consistent within uncertainties. The quark PDFs are essentially unaffected 
by the inclusion of the LHCb charm data and are thus not shown here. 

Since the resulting PDF set from the inclusion of the LHCb data into NNPDF3.0 has been 
unweighted to a a LHAPDF6 grid, it can be easily used both for the predictions of heavy quark 
production at 13 TeV at LHCb, presented in Sect. and for the prompt neutrino cross-sections 
relevant for IceCube in Sect [H 

It is interesting to assess how the results of this analysis compare to those of the PRDSA 
study 37 . Note that the two analysis use rather different methodologies (HERA-only fit versus 


global fit, HERAfitter versus NNPDE reweighting), and given that this is the first time that 
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Figure 12: Left: The NNPDF3.0 NLO small-a: gluon, evaluated at Q = 2 GeV, comparing the global 
fit result with with the new gluon obtained from the inclusion of the LHCb charm production data. In 
the latter case, we show both the reweighted {rwg) and the unweighted [unw) results. Right: comparison 
of percentage PDF uncertainties for the NNPDF3.0 gluon with and without the inclusion of the LHCb 
data, computed also at Q = 2 GeV, that illustrate the reduction of PDF uncertainties for x < 10“^. 


forward charm data is used in a PDF fit, it is important assess the robustness of the results by 
performing a cross-check. Since the PROSA analysis is performed in the FFN nj = 3 scheme, we 
have constructed a FFN n/ = 3 version of the NNPDF3.0-|-LHCb NLO set using APFEL j^. The 
results of this comparison are shown in Fig. 13, where we show the gluon PDF at = 10 GeV^ 


in the FFN scheme with Nf = 3, In the PROSA case, we show the results both in the HERA-only 
fit and in the HERA-l-LHCb fitj^ The lower panel compares the relative PDF uncertainties 
in each case. As can be seen, there is good agreement both between central values (the two 
gluons agree within their one-sigma band) and especially between PDF uncertainties, which is 
a non-trivial verification of the two analyses. 

Finally, let us compare the resulting gluon PDF in this analysis with those of other recent 
PDF fits. In Fig. 14 we compare the NNPDF3.0-|-LHCb gluon PDF at = 4 GeV^ with the 


GT14 and MMHT14 results (left plot), and to the ABM12 84 and HERAPDF2.0 results 
(right plot). In the case of HERAPDF2.0, both the experimental, model and parametrization 
uncertainties are included. In the case of ABM12, the n/ = 4 set has been adopted. From 
Fig. 14 we note that the NNPDF3.0-|-LHCb central value is close to the CT14 result, but with 
much smaller uncertainties, while the MMHT14 gluon is substantially larger at small-x. From 
the comparison with ABM12 we find reasonable agreement for x < 10“^, while HERAPDF2.0 
predicts a much smaller (negative gluon), though consistent with the NNPDF3.0-|-LHGb result 
within the PDF large uncertainties. 


4 Predictions for 13 TeV and for the 13/7 TeV ratio 

In this section we provide predictions for D and B production within the LHGb acceptance 
at 13 TeV. We also provide predictions for the ratio of differential cross-sections between 13 
and 7 TeV. Our predictions are have been computed using the POWHEG and aMC@NLO 
calculations with the improved NNPDF3.0-|-LHGb PDF set constructed in Sect.j^ and can be 
used to compare with the upcoming Run II measurements at LHCb. Using the theoretical value 
of the ratio between inclusive fiducial cross-sections at 13 and 7 TeV, and the LHCb 7 TeV data 

^We thank Katerina Lipka for providing us this plot, which compares the PROSA and NNPDF results. 
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Figure 13: The gluon PDF at = 10 GeV^ in the FFN scheme with Nf = 3, comparing the results of 
this work with those of the PROSA analysis. In the latter case, we show the results both in the HERA-only 
fit and in the HERA+LHCb fit. The lower panel compares the relative PDF uncertainties in each case. 


{R 13 / 7 ), we also provide predictions for B and D mesons in fiducial cross-sections at 13 TeV. A 
tabulation of our results is provided in Appendix [A] and predictions for different binning choices 
and other meson species are available from the authors on request]^ 


4.1 Forward heavy quark production at 13 TeV 

First of all, we provide theory predictions required to compare with the upcoming LHCb data 
on charm and bottom production which will be collected at 13 TeV. Our results are presented 
according to the binning scheme adopted in the 7 TeV measurements 35,36 , with the exception 


that a slightly finer binning for the charm predictions is chosen at low pT and the high pT range 
is slightly extended. For all predictions, the uncertainty due to scales, PDFs, and the heavy 
quark mass is provided as a sum in quadrature. 

the double differential distributions for mesons at 13 TeV are shown for 


In Fig. 15 


both a central and a forward rapidity bin within the LHCb acceptance. The central value and 
total uncertainty of both POWHEG and aMC@NLO calculations are provided. This comparison 
demonstrates that there is good agreement between the two calculations, both in terms of central 
values and in terms of the total uncertainty band — agreement also holds for other D mesons 
and rapidity regions, which are not shown here. Thanks to using the improved NNPDF3.0 PDFs 
with 7 TeV LHCb data, PDF uncertainties turn out to be moderate even at 13 TeV, with scale 
variations being the dominant source of theoretical uncertainty. 


®Very recently, the LHCb 13 TeV charm production measurements have been presented 


87 


The LHCb 


publication includes a detailed comparison between data and the theoretical predictions presented in this work, 
showing good agreement within uncertainties. This agreement for the 13 TeV data provides further validation of 
the robustness of our approach. 
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Figure 14: The NNPDF3.0+LHCb gluon PDF at = 4 QeV^ compared with CT14 and MMHT14 
(left plot), and to ABM12 and HERAPDF2.0 (right plot). In the case of HERAPDF2.0, both the 
experimental, model and parametrization uncertainties are included. 



Figure 15: The double-differential distribution, (Pa{D)/dydpT, for the production of mesons at LHCb 
for a centre-of-mass energy of 13 TeV. We show representative results for the central (2.0 < y < 2.5) and 
forward (3.5 < y < 4.0) regions. We compare the POWHEG and aMC@NLO calculations, using the 
NNPDF3.0-|-LHCb PDF set. For both calculations, the theory uncertainty band is computed adding in 
quadrature scales, PDF and charm mass uncertainties. 


The corresponding comparison for mesons is shown in Fig. 16 As in the case of the 


charm, there is excellent agreement between the POWHEG and aMG@NLO calculations within 
the LHCb acceptance. 

The tabulation of the results shown in Figs. 15 and are provided in Appendix in 
particular in Tables]^ (for mesons) and(for mesons). 


4.2 Predictions for the ratio between the 13 and 7 TeV cross-sections 

In addition to differential cross section measurements, it will also become possible to measure the 
ratio of differential cross sections performed at 13 and 7 TeV when the 13 TeV data is available. 
As discussed in Ref. , measurements of the ratio of cross-sections at different centre-of-mass 
energies are well motivated as many theoretical uncertainties, such as scale uncertainties, mass 
dependence, and fragmentation/branching fractions cancel in the ratio to a good approximation. 
In addition, many experimental uncertainties also cancel in such ratios which allows stringent 


21 






















NNPDF3.0+LHCb, 6° mesons, 2.0 < y < 2.5 


NNPDF3.0+LHCb, B“ mesons, 3.5 < y < 4.0 




Figure 16: Same as Fig. 15 for mesons. 


tests of the Standard Model to be performed. The relevance of the ratio of 13 over 7 TeV heavy 
quark production cross-sections at LHCb for PDF studies has also been recently emphasised in 
Ref. [89| , in a study of the various theoretical uncertainties associated to charm and bottom 
production in the forward region. 

On the other hand, PDF uncertainties do not cancel completely, because of the different 
kinematical range covered by the measurements at the two centre-of-mass energies, and thus 
these ratio measurements provide in principle useful PDF discrimination power. This idea has 
been implemented already by a number of LHC analyses, like the ATLAS measurement of the 
ratio of 7 TeV over 2.76 TeV jet cross-sections 83 and the CMS measurement of the ratio of 8 


TeV over 7 TeV Drell-Yan distributions 88 


In Fig. 17 we show the predictions for the ratio of differential cross-sections for production 


between 13 TeV and 7 TeV, defined as 


R 


13/7 


D D^ cf a13 TeV) 
- dM - 


d^a{D^){y^,p^,7TeV) 

dy^dpS 


( 4 ) 


where the same binning as in the 7 TeV LHCb measurement has been assumed. In the left plot 
we show the results computed with POWHEG and the NNPDF3.0-|-LHCb PDF set, for each of 
the bins of the 7 TeV measurement (data points are ordered in increasing bins of rapidity, and 
within each of these five rapidity bins, in increasing bins of pt)- The central value of the ratio 
'^ 13/7 between 1.20 and 2.2 for increasing values of px and more forward rapidity bins, 

where the opening of phase space between 13 TeV with respect to 7 TeV is more important. 


In the left plot of Fig. 17 we have separated the total theory uncertainty into the individual 


contributions from scales, PDFs and charm mass to highlight their importance. We see that 
the total uncertainty in y varies between 10% and 30%, depending on the specific bin, and 

that scale variation is found to dominate the total uncertainty in Note however the 

substantial cancellation of scale uncertainties as compared to the absolute differential cross- 
sections shown in Fig. 15 In Appendix we provide a tabulation of the results of Fig. 0 


which will be useful for comparison if the ratio y is measured in the upcoming LHCb 13 TeV 

analysis. In the same appendix we also quantify the reduction of PDF uncertainties in R^y^ 
by comparing the predictions using the original NNPDF3.0 set with our baseline predictions 
obtained with NNPDF3.0-|-LHCb. The substantial reduction of PDF uncertainties in R 


thanks to the constraints from the 7 TeV normalised charm cross-sections derived in Sect. 


D° 

13/7’ 
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Figure 17: Left plot: predictions for the ratio of differential cross-sections Eq. for the pro¬ 

duction of mesons between 13 TeV and 7 TeV, computed using POWHEG and NNPDF3.0-|-LHCb. 
Results are ordered in increasing bins in rapidity, in within each, in increasing bins of px- The total theo¬ 
retical uncertainty in the ratio is decomposed into its various sources: scale, PDF and charm quark mass 
variations. Right plot: comparison of the predictions for Rf ^/7 between POWHEG and aMC@NLO, 
for central values and for the total theory uncertainties. 


improve the robustness of our theory prediction for 7 - Conversely, the measnrement of 7 
should provide important PDF discrimination power, and it would be interesting to verify the 
consistency of the constraints on the small-x gluon from 7 from those that we have derived 
from the normalised 7 TeV data. 

To validate the cancellation of the theoretical systematics in the POWHEG calculation, we 
have also computed the ratio with the aMG@NLO calculation. The comparison of these two 


calculations, including their total uncertainties, is shown in the right plot of Fig. 17 Reasonable 
agreement is found, both for the central values and for the uncertainties. In particnlar, for most 
of the bins, the central predictions for R (^°7 agree within 10% at most. This agreement should 
be considered satisfactory especially taking into account the very large theory uncertainties in 
the absolute distributions. 

Next we provide the corresponding predictions for the ratio of B meson differential distribu¬ 
tions between 13 TeV and 7 TeV, dehned as 




B TeV) (fcj{B^){yB,p^,lT(N) 


dy^dp^ 


dy^dp^ 


(5) 


for the case of B^ mesons, which we choose for illnstrative pnrposes. In Fig. 18 we show the 
theoretical predictions for the ratio R ^/7 computed with POWHEG using NNPDF3.0-|-LHCb 
for two representative bins in rapidity, one central (left plot) and one forward (right plot), as 
a function of p^. The total theory uncertainty (hatched band) is compared with the scale 
uncertainty (solid band). We have verified that the results for R ^'/7 obtained with aMG@NLO 
are fnlly consistent the POWHEG calculation. In the results of Fig. 18, the same binning as in 
the 7 TeV measurement has been used 36 


From Fig. 


18 


we see that R ^/7 varies between 1.3 at central rapidities at low pp to almost 


at forward rapidities and large pp, for the same reasons as The total uncertainty in 


Rf^^j ranges between 5 and 10%, depending on the specific bin, and is dominated by the scale 
nncertainty (but only due to using the improved NNPDF3.0-|-LHCb set). As in the case of 
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NNPDF3.0+LHCb, B“ mesons, 2.0 < y < 2.5 NNPDF3.0+LHCb, B“ mesons, 3.5 < y < 4.0 




Figure 18: Theoretical predictions for the ratio Eq. |5| between meson distributions between 

13 and 7 TeV. Results have been computed with POWHETG using NNPDF3.0+LHCb. We show the 
predictions for two representative bins in rapidity, one central (left plot) and the other forward (right 
plot), as a function of p!^. The total theory uncertainty (hatched band) is compared with the scale 
uncertainty (solid band). 


charm production, in Appendix A 
compare the the upcoming LHCB" measurement, 
meson species are available upon request. 


we tabulate our predictions for that can be used to 

The corresponding predictions for other B 


4.3 Predictions for inclusive fiducial cross-sections at 13 TeV 

In addition to the double differential distributions, it is also useful to provide predictions for the 
charm and bottom inclusive cross-section, that is, the cross-sections measured within the full 
LHCb fiducial region. In the case of D mesons, the fiducial region is defined as 

0<p^<8GeV, 2.0<y^<4.5, (6) 

while the corresponding fiducial region for the production of B mesons is defined by 

0 < pf < 40 GeV, 2.0 < < 4.5. (7) 


In order to compute the 13 TeV predictions for the charm and bottom inclusive cross-sections 
in the fiducial region, there are two possible strategies that can be adopted, namely 


integrating the POWHEG calculation for the absolute double differential cross-sections, 
shown in Fig. 15, for the acceptance in Eq. (|^, or instead 

using the theoretical predictions for the ratios ^^*4 Ri^/j to rescale the corresponding 


7 TeV LHCb inclusive measurements reported in 35,36 


The main advantage of the second option is that theoretical uncertainties are substantially 
reduced in the ratios R^^^j and Rf^^^ compared to the absolute cross-sections, allowing a rea¬ 
sonably accurate extrapolation for the 13 TeV inclusive cross-sections, with precision comparable 
to that expected for the corresponding experimental measurement. 

Let us illustrate how the two strategies compare in the case of D meson production. For 
simplicity, we will show the results for mesons but the same ideas apply to the other D 
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mesons. In this case, the prediction for the inclusive ratio, with the total associated theory 
uncertainty, is given by 

-^13/7_o.29(20.5%) ' 


This can be combined with the 7 TeV LHCb inclusive measurement 36 in the fiducial region 
for mesons 

o-|^ev(LHCb,incl) = 1661 ± 129 (±7.8%) /xb , (9) 

to obtain an accurate prediction for the corresponding 13 TeV inclusive cross-section in the same 
fiducial region. This leads to 


^^f 3 Tev(th,incl) = (T^ev(LHCb,incl) • 72 ^ 3 / 7 (th, incl) = 2236(^521(23%) ’ 






,-1-308 (14%) 


( 10 ) 


where the theoretical uncertainty from is slightly larger than that of the 7 TeV mea¬ 
surement, and dominates the precision of the prediction for performed in this way. In 

Eq. (10) we have added in quadrature the theory uncertainties from 7 with the experimental 
uncertainties of the LHCb measurement. 

In Table the prediction for the inclusive cross-section obtained using the 7 TeV 

measurement and the calculation of 72^^^ is compared to the corresponding result computed 
from the integral of the absolute differential distributions. The advantage of the ratio strategy is 
apparent; when integrating the absolute distributions, the prediction is affected by large theory 
uncertainties up to 200 % which render the comparison with the much more accurate experimen¬ 
tal measurement not very informative. On the other hand, our prediction obtained using 72^^.^, 
has a 10-20% accuracy, comparable to that of the upcoming Run II LHCb measurement, and 
therefore should provide interesting information for the comparison between data and theory 
in a hitherto unexplored kinematical region. In Table we also provide the predictions for the 
inclusive charm pair production cross-section using the two methods, obtained from rescaling 
the meson-level result by the branching fraction of charm into mesons, 

0'i3Tev(th, incl) = o-(^Tev(th, mcl)/ ( 2 / (c ^ 77°)) . (11) 

This prediction is useful to compare with parton-level predictions of charm production, which 
do not account for the fragmentation of charm quarks into D mesons. 


13 TeV 

77° 

CC 

<7'i3Tev(th, incl)(/x 6 ) (from ratio) 

5 , 5 , 0 ^ +308(14%) 
ZZOD 52i(23%) 

1970 +249(13%) 
ly'y -447(23%) 

<7'i3Tev(th, incl)(^ 6 ) (from abs) 

+2082(190%) 

-896(82%) 

+1843(190%) 

-793(82%) 


Table 2: Predictions for the inclusive 77° production cross-section in the hducial region Eq. (|^ at 13 
TeV using the two methods discussed in the text (integrating the absolute distributions and rescaling the 
7 TeV LHCb measurement with the ratio Predictions are also provided for the corresponding cc 

cross-sections using Eq. [m 

The same strategies can be applied to obtain accurate predictions for the inclusive B meson 
production cross-sections at 13 TeV in the fiducial region defined by Eq. ([^. For simplicity we 
restrict ourselves to 77° mesons, though the same method also applies to all other B mesons that 
will be measured at Run H. The first method, integrating the absolute differential cross-sections 


from Fig. 16 in this fiducial region leads to the following prediction 

< 7 f 3 TeV(til, incl)(/x 6 )(from abs) = 55.07 ^2a76 (3L7%) ■ 


( 12 ) 
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Now, using the prediction for the ratio of inclusive cross-sections between 13 and 7 TeV for 
mesons, 

Rfs/jith, incl) = 1.84 > (13) 

to rescale the 7 TeV LHCb measurements in this fiducial region, 

CT’^gV(LHCb) ™cl) = 38.1 ± 6.0 (±15.6%) /rb , (14) 

we obtain the following prediction for the 13 TeV hducial production cross-section 

0'i3Tev(th5 incl)(from rat) = iT^gv(LllCl’) ™cl) • iif 3 '’^ 7 (th, incl) = 70.02 03 (i 7 2%) ' (1^) 


In the above procedure, the theoretical uncertainties from 7?^^^ and the experimental uncertain¬ 
ties from the 7 TeV measurement have been added in quadrature. In this case, the advantage of 
using more marked: as the theoretical uncertainties of the ratio are smaller than 

those of the 7 TeV LHCb inclusive measurement, the extrapolation from 7 to 13 TeV is essen¬ 
tially limited by the precision of the 7 TeV cross-section, with very small theoretical uncertainty 
in the procedure. Note that using the ratio strategy our theoretical prediction for uisTeV leads to 
a prediction with uncertainties which are around three times smaller as compared to the predic¬ 
tion obtained from the integration of the absolute distributions, Eq. (12). Similar improvements 
can be observed for other meson species. Note also that in the case of B mesons, the fragmen¬ 
tation is essentially the same for all the meson types, and thus the same rescaling Eq. (13) can 
be applied to all the B meson species. For example, we find i?^^y(th, incl) = 77^3°^.^(th, incl) to 
the precision provided in Eq. (13). 

In summary, in this section we have provided accurate predictions for the 13 TeV fiducial 
cross-sections for the production of D and B mesons at LHCb, using the ratios R^^^j -^ 13/7 

to extrapolate the 7 TeV measurements. The robustness of this extrapolation is illustrated by 
the fact that, upon rescaling by the ratio, the corresponding 13 TeV prediction has uncertainties 
which are at most two times larger than than the precision of the 7 TeV data. Note that the 
predictions from the absolute distributions have significantly larger uncertainties as compared 
to the foreseen prediction of the 13 TeV uncertainties, particularly in the case of charm, where 
theory uncertainties for the fiducial cross-section can be as large as ~ 200% (see Table . 


5 QCD predictions for charm-induced neutrino production 

The dominant background for the detection of ultra-high-energy neutrinos from astrophysical 
sources in experiments like IceCube arises from the flux of neutrinos originating from the prompt 
decay of energetic charmed mesons produced in cosmic ray collisions in the upper atmosphere. 
We now provide state-of-the-art pQCD predictions for the cross-sections of charm-induced neu¬ 
trino production. These cross-sections are an important ingredient of the full calculation of 
prompt neutrino event rates at IceCube, which is beyond the scope of this paper. 

As compared to previous works [^13|, here we want to fully exploit the flexibility of our 
approach for the computation of the charm production cross-sections, based on NLO Monte 
Carlo event generators. We can derive a robust prediction for the primary neutrino flux arising 
from the decays of charmed mesons produced in cosmic ray collisions from pQCD, eliminating the 
need of model assumptions, and being able to estimate all the associated sources of theoretical 
uncertainties in our calculation. Being fully differential, our calculation of the prompt neutrino 
flux can be processed in cascade codes and in neutrino telescopes detector simulation software 
with arbitrary selection cuts. 
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To achieve this goal, using the results of Sects. and we have computed 

da{pN —> i/X; E] E^) 


( 16 ) 


that is, the differential cross-section for the production of neutrinos from the decays of charmed 
hadrons in proton-nucleon collisions, as a function of the neutrino energy Ey, for different values 
of the incoming cosmic ray energy Ef\ 


To compute the neutrino energy distribution, Eq. (16), using POWHEG and aMC@NLO 


charmed hadrons are first decayed using the PythiaS shower, summing over all hadron species 
and neutrino flavours. Subsequently, a Lorentz boost is applied for the conversion of the neutrino 
energy distribution from the centre-of-mass frame, where the prediction of MC event generators 
is provided, to the laboratory frame. The magnitude of this boost is determined by the incoming 
cosmic ray energy. 

Results have been computed for a number of values of the incoming cosmic ray energy E 
between E = 10^ GeV to E = 100 PeV, corresponding to centre of mass energies \/s = 
ranging from 44 GeV to 14 TeV. As discussed before, we emphasize the overlap between the 
kinematic region crucial for neutrino telescopes and that of the LHGb charm production data. 

The fact that cosmic rays collide with air nucleus rather than with isolated (isoscalar) nucle¬ 
ons can be accounted for by rescaling the cross-section for pN collisions with the mean atomic 
number of air nuclei (A) ~ 14.5, that is, to good approximation we can write 


a{pA —)> ccX) ~ (A) • a{pN —)> ccX). 


(17) 


Eq. ( |17[ ) assumes that nuclei can be treated as an incoherent sum of their protons and neutrons, 
and that nuclear corrections to the nucleon PDFs can be neglected as compared other theoretical 
uncertainties in the calculation. 

The assumption of neglecting nuclear shadowing in charm production is justified by the recent 
GMS measurements of B mesons in proton-lead collisions at ^/snn = 5 TeV 90 , which cover 
a similar kinematical range as for charm production in cosmic rays, and that show no evidence 
for suppression induced by nuclear PDFs. Moreover, available sets of nuclear PDFs 91-93 


are unconstrained at small-x due to the absence of experimental data, and thus cannot be 
used reliably in our calculation. In addition, a recent calculation of forward D production at 
= 5 TeV incorporating the EPS09 nuclear PDF modifications 71 indicates that a cross 
section suppression of at most ~ 10% can expected in proton-lead collisions, within substantial 
uncertainties. 

In our approach the production of charm quarks, their hadronisation into charmed mesons 
and their subsequent decays into neutrinos are completely accounted for in the matrix element 
calculation matched to the parton shower. We can therefore obtain exact results for the various 
differential distributions relevant for prompt neutrino production. We emphasise that the the 
modelling of charm production and decay in PythiaS has been validated by LEP data as well 


as hadron collider data, see Refs. 166,94 and references therein. 


These differential cross-sections Eq. (16) have been computed in a range of values of E 
and Ey and then suitably interpolated. For each point in {E^E^,), we have determined the 
relevant theoretical uncertainties from scales, PDFs, and rric variations. Our calculations use 
the improved NNPDF3.0-|-LHGb which includes the constraints from the 7 TeV charm data. A 


°Eq. (161 accounts only for the flux of primary prompt neutrinos, those produced in the first interaction of the 
cosmic ray with air nuclei. To compute the complete flux one should also include the contribution from secondary 
production solving the cascade equations. 
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Figure 19: The differential cross-section for the production of neutrinos from charm decay in pp collisions, 
Eq.[^ as a function of the neutrino energy, computed with POWHEG. The results are provided for two 
values of the incoming cosmic ray energy, E = 10^ GeV (left plot) and E = 10® GeV (right plot). The 
input PDF set is NNPDF3.0NLO-|-LHCb. We show the central prediction as well as the scale, PDF and 
TOc uncertainties, as well as the overall theoretical uncertainty band computed from adding in quadrature 
the three independent theory errors. 


representative sample of our predictions are provided in Fig. 19 where we show the differential 


cross-section for the production of neutrinos from charm decay in pp collisions, Eq. (16), as a 


function of the neutrino energy, computed with the POWHEG calculation. Results are shown 
for two values of the incoming cosmic ray energy, E = 10^ GeV (left plot) and E = 10® GeV 
(right plot). We show the central prediction as well as the individual contributions from scale, 
PDF and nic uncertainties, as well as the overall theoretical uncertainty band computed from 
adding these uncertainties in quadrature. We see that at the highest energies, E = 10® GeV, 
the total uncertainty band is dominated by scale variations, while PDF uncertainties are under 
control thanks to the constraints from the LHCb charm production data. We stress that while 
NLO QCD scale uncertainties are still large, up to a factor three, recent work towards the NNLO 
differential distributions for heavy quark production will provide a reduction of these 

higher-order uncertainties. 


A powerful cross-check of the robustness of the predictions shown in Fig. 19 is provided 
by the fact that comparable results are obtained using either POWHEG or aMG@NLO, both 
for the central prediction and for the upper and lower ranges of the total theory uncertainty 


band, as shown in Fig. 20, when the same theory settings are used in the two calculations. 


Let us emphasise that two completely independent codes are used, with different underlying 
matrix element calculations and different matching to the parton showers, so this agreement 
is an indication of the robustness of the pQCD predictions for the charm-induced neutrino 
production cross-sections presented here. 

It is also interesting to study the dependence of our results for the charm-induced neutrino 
production cross-sections as a function of the incoming cosmic ray energy E. In Fig. 21 


we 


represent the differential cross-section for neutrino production in charm decays, Eq. (16), for 
different values of E, as a function of the ratio between the neutrino energy E^^ and the cosmic 
ray energy, z = E^/Ep, that is. 


da{pN —)• vX] E; Ey = zE) 


Ey 


z = 


(18) 


dz ’ ~ E ’ 

which allows to compare the increase of the neutrino production cross-section, due to the larger 
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Figure 20: Same as the right plot of Fig. 19 now comparing the predictions of POWHEG with those 
of aMC@NLO. The same theory settings are used in the two calculations. Only the central curve total 
theory uncertainty bands are shown for the predictions obtained with the two event generators. 


value of E, for the same value of z, the ratio of the neutrino energy over the incoming cosmic 
ray energy. In Fig. 21 results are shown for E = 10^ and E = 10® GeV (both central values 


and total theoretical uncertainty) and then for E = 10® and E = 10® GeV (only central values). 
Note how the cross-sections fall steeply as one approaches the kinematical boundary, z —>■ 1. 

Note that in pQGD, the correct expression for representing the dependence of E of the 
prompt neutrino production cross-section is given by Eq. (18), shown in Fig. 21 Previous 
works, for example |13| , present their calculations of the charm production cross-section as 
da/dxc, where Xc = Ec/E, the ratio of produced charm quark energy over the incoming proton 
energy. However, the charm quark energy is only well dehned at leading order, beyond which 
this is not true, and moreover is not accessible experimentally. Therefore, a robust comparison 
of theoretical calculations should always be presented at the level of the physical D production 
cross-section. Alternatively, one might rescale by the charm branching fraction as in Eq. 
but this approximation is only valid for relatively inclusive observables. 


Einally, let us mention that our calculations for Eq. (16), illustrated in Eigs. 19 and|20[ are 
available for a wide range of E and Ey values in the format of interpolated tables that can be 
used as input for calculations of the prompt neutrino flux at IceCube, and are available from 
the authors upon request. 


6 Summary and outlook 

In this work we have performed a detailed study of charm and bottom production in the forward 
region, based on state-of-the-art pQCD with NLO calculations matched to parton showers. Our 
motivation was to provide a robust estimate of the theoretical uncertainties associated to the 
prompt neutrino flux at neutrino telescopes like IceCube, which is the dominant background for 
the detection of astrophysical neutrinos. 
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Figure 21: The dependence on the incoming cosmic ray energy E of the prompt neutrino production 
cross-section da/dE^, plotted as a function of z = E^/E, Eq. (18), which allows to compare calculations 
for different values of E. Results are shown for E = 10^ and E = 10® GeV (both central values and total 
theoretical uncertainty) and then for E = 10® and E = 10® GeV (only central values). The cross-sections 
fall steeply as one approaches the kinematical boundary z —>■ 1. 


Our strategy was based on the careful validation of the pQCD calculations with the LHCb 
charm and bottom production data at 7 TeV, which cover the same kinematical region as that 
relevant for the production of prompt neutrinos at IceCube. We found that, with a suitable 
normalisation of the differential distributions, it is possible to include the 7 TeV D meson data 
from LHCb in order to significantly constrain the poorly known small-x gluon. Being able to 
include the LHCb charm measurements in a global NLO PDF fit further enhances our confidence 
of the applicability of pQCD to provide predictions for the prompt neutrino flux. These improved 
PDFs, NNPDF3.0-|-LHCb, which include the information from the LHCb charm data, are then 
used to construct the predictions for charm and bottom production at LHCb for the recently 
started Run H with a centre-of-mass energy of 13 TeV, as well as for the ratio of 13 over 7 TeV 
cross-sections. 

Our main result for the cross-sections of the production of prompt neutrinos in charmed me¬ 


son decays originating from cosmic ray collisions in the atmosphere is summarised in Figs. 19 20 


and 21 The main difference as compared to previous calculations of the prompt neutrino flux 


is that our approach has been fully validated with the recent LHCb differential measurements 
on charm production, and that the input PDF set used in our calculations is one that already 
includes the constraints from the LHCb charm data. We would like to emphasise that our cal¬ 
culations, both for the central values and for the various theoretical uncertainties, have been 
carefully benchmarked using three independent codes. 

The main results of our study can be summarised as follows: 


• pQCD predictions for charm and bottom production in the forward region are consistent 
with the recent LHCb 7 TeV measurements within theoretical uncertainties. Predictions 
obtained with three different codes, two Monte Carlo parton shower programs, aMC@NLO 
and POWHEG, and one semi-analytical calculation, FONLL, yield comparable results. 
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both for the central value and for the total uncertainty. 

• It is possible to include the LHCb charm data in the NNPDF3.0 NLO global analysis, 
achieving a substantial reduction of the PDF uncertainties on the poorly known small-x 
gluon. In order to reduce the large scale uncertainties of the NLO calculation, the LHCb 
data have been normalised to a hxed reference bin. 


• Run II of the LHC has just started, and the LHCb experiment will soon measure charm and 
bottom production in the forward region at 13 TeV, which will further explore the low-x 
region of gluon PDF providing unique information on the structure of the proton. We have 
thus provided predictions for charm and bottom production at 13 TeV, as well as for the 
ratio of differential cross-sections between 13 and 7 TeV. These new measurements, both 
the 13 TeV (normalised) differential distributions and the 13 over 7 TeV cross-section ratio, 
offer new possibilities for PDF constraints, in particular thanks to the extended coverage 
at small-a: as compared to the 7 TeV measurements. 

• Using the theory prediction for the ratio of inclusive fiducial cross-sections R 13/7 combined 
with the corresponding LHCb 7 TeV measurements, we are able to provide a prediction for 
the 13 TeV fiducial cross-section with substantially reduced uncertainties as that compared 
to the prediction from the NLO QCD calculation. 


We have provided QCD predictions for the differential cross-sections for the production of 
neutrinos from charm decay in pA collisions, Eq. (16), using two independent NLO Monte 
Carlo generators, across a wide range of incoming cosmic ray energies E, and accounting 
for all relevant theory uncertainties. 


It will be interesting to compare the upcoming 13 TeV LHCb measurements with the pre¬ 
dictions presented in this paper. In particular, one should verify that the constraints on the 
small-x gluon obtained from the inclusion on the PDF fit of the measurement of the ratio R 13/7 
of differential distributions are consistent with those that have been obtained from the 7 TeV 
normalised charm production cross-sections. Likewise, comparing the inclusive fiducial cross- 
sections at 13 TeV with our predictions based on R 13/7 will be an important test of the validity 
of QCD calculations in this new kinematical region. 

It is beyond the scope of this paper to explore quantitatively the implications of our calcu¬ 
lations for the recent IceCube measurements of ultra high energy neutrinos. Our results for the 
charm-induced neutrino cross-sections as a function of E and Ei, are available in the form of 
interpolated grids. This information can be used as input in a full calculation to derive robust 
predictions for the rates of prompt neutrino events expected at IceCube. 
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A Predictions for charm and bottom production at 13 TeV 


In this appendix we provide a tabulation of our predictions for charm and bottom production 
in LHCb at TeV, presented in Sect. |4.l] as well as for the ratio of cross-sections between 13 TeV 
and 7 TeV, discussed in Sect. 4.2 For simplicity, we will restrict ourselves to the POWHEG 


results, since we have established from the comparison with aMC@NLO in Figs. 15 and |16| that 
the two calculations yield similar results. 


A.l Predictions for differential distributions at 13 TeV 

First of all, in Table we provide the predictions for the differential cross-sections for 


production at 13 TeV in the LHCb acceptance, corresponding to the results in Fig. 15 These 


results have been obtained with POWHEG using NNPDF3.0-|-LHCb as input PDF. For each 
bin, we provide the central value and the total theoretical uncertainty. To take into account 
the increased statistics that the Run H measurement will benefit from, we have used in this 
tabulation an optimised binning as compared to the 7 TeV results. First of all, we have used a 
hner binning at low pT (the region which is most sensitive to the gluon PDF) and extended our 
predictions up top^ = 30 GeV (where theoretical uncertainties are smallest). The corresponding 
predictions for any other choice of binning in as well as for the other D meson species 

are available from the authors upon request. As in the case of the 7 TeV results, absolute cross- 
sections are affected by substantial theoretical uncertainties, in particular due to the large scale 
variations of the NLO computation. On the other hand, PDF uncertainties are now subdominant 
for all values of and p!^, thanks for the constraints from the 7 TeV LHCb charm data. 

The corresponding predictions for the differential cross-sections of the production of 
mesons at LHCb Run H are shown in Table which is the analog of Table for charm pro¬ 
duction. These predictions were represented graphically (and compared to the aMC@NLO 


calculation) in Fig. 16 In this case we have assumed the same binning as in the 7 TeV measure¬ 


ment. As compared to the 13 TeV charm predictions, the higher scales and the larger values of 
Bjorken-x probed in the case of bottom production result in reduced theory uncertainties. The 
total uncertainty is around 50%, with differences depending on the specific bin, and is again 
dominated by scale uncertainties. 


A.2 Predictions for the ratio R 


13/7 


Next we turn to the predictions for the ratio R 13/7 of the differential distributions for heavy 
quark production at LHCb between 13 TeV and 7 TeV, Eqns. Q and ([^, discussed in Sect. 4.2 
In Table we show the ratio R 13/7 for mesons using POWHEG and NNPDF3.0-|-LHCb. 


These prmictions were represented graphically in Fig.j^ We provide the central value of R 13/7 
and the total theoretical uncertainty, which as can be seen from Fig. HZl arises predominantly 
due to scale variations. When evaluating Eq. Q, scale variations, charm mass variations and 
PDF variations are considered to be fully correlated between 13 and 7 TeV. 

In order to evaluate the impact of the reduction of PDF uncertainties on the observable Rf^/ 7 , 
that has been achieved by including in NNPDF3.0 the LHCb 7 TeV charm production data, it is 
useful to compare with the corresponding predictions with the original NNPDF3.0 set. With this 
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(Mfe/GeV) 


P? (GeV) 



2.0 

-2.5 

2.5 

-3.0 

3.0 

-3.5 

3.5 

-4.0 

4.0- 

4.5 

0.0-0.25 

16.7 

+ 41.4 

- 15.8 

16.4 

+ 40.6 

- 15.3 

16.2 

+ 41.0 

- 14.9 

15.0 

+ 38.0 

- 13.8 

13.5 

+ 34.0 

- 12.2 

0.25-0.5 

45.4 

+ 109.8 
- 43.0 

45.9 

+ 114.6 
- 42.8 

43.6 

+ 108.7 
- 40.2 

41.5 

+ 103.5 
- 37.6 

38.6 

+ 96.4 

- 34.1 

0.5-0.75 

66.0 

+ 156.6 
- 61.4 

65.2 

+ 157.6 
- 60.0 

63.0 

+ 153.0 
- 57.1 

59.9 

+ 146.3 

- 53.8 

55.4 

+ 136.1 
- 48.8 

0.75- 1.0 

77.3 

+ 178.6 
- 71.7 

75.1 

+ 174.6 
- 68.5 

73.5 

+ 171.5 
- 65.7 

70.5 

+ 167.5 

- 62.8 

63.0 

+ 148.6 
- 55.2 

1.0- 1.25 

82.2 

+ 181.7 
- 75.0 

80.1 

+ 178.5 
- 72.2 

76.1 

+ 171.0 
- 67.5 

72.4 

+ 164.2 
- 63.4 

66.8 

+ 152.9 
- 57.9 

1.25 - 1.5 

79.6 

+ 168.0 
- 71.8 

77.7 

+ 165.2 
- 69.2 

73.8 

+ 155.5 
- 64.8 

69.9 

+ 150.3 
- 60.2 

63.7 

+ 138.2 
- 54.6 

1.5- 1.75 

76.2 

+ 152.2 
- 67.7 

73.4 

+ 148.8 
- 64.1 

69.9 

+ 143.0 
- 60.0 

65.8 

+ 137.1 

- 56.0 

59.3 

+ 124.0 
- 49.8 

1.75-2.0 

69.7 

+ 132.7 
- 60.9 

67.7 

+ 129.5 
- 58.5 

64.8 

+ 127.6 
- 54.8 

59.7 

+ 117.2 

- 49.2 

53.1 

+ 106.5 
- 44.0 

2.0-2.25 

62.7 

+ 112.3 
- 53.8 

59.8 

+ 108.0 
- 50.1 

57.5 

+ 107.0 
- 47.4 

52.0 

+ 97.1 

- 42.6 

46.9 

+ 89.1 

- 37.6 

2.25 - 2.5 

55.3 

+ 95.1 

- 46.3 

53.1 

+ 91.5 

- 43.8 

50.2 

+ 88.4 

- 40.8 

46.3 

+ 82.2 

- 36.8 

40.2 

+ 72.6 

- 31.8 

2.5-2.75 

49.7 

+ 82.8 

- 41.1 

46.5 

+ 77.1 

- 37.7 

43.7 

+ 73.2 

- 34.7 

39.3 

+ 67.2 

- 30.8 

34.7 

+ 58.8 

- 26.7 

2.75-3.0 

43.4 

+ 68.7 

- 35.0 

41.2 

+ 66.2 

- 32.5 

37.9 

+ 60.6 

- 29.2 

34.0 

+ 56.2 

- 25.8 

29.3 

+ 49.2 

- 21.8 

3.0-3.5 

35.7 

+ 53.6 

- 28.1 

33.5 

+ 50.6 

- 25.7 

30.7 

+ 46.8 

- 23.1 

27.6 

+ 42.6 

- 20.4 

22.9 

+ 35.9 

- 16.8 

3.5-4.0 

26.9 

+ 37.5 

- 20.2 

25.0 

+ 34.8 

- 18.3 

22.9 

+ 32.4 

- 16.5 

20.0 

+ 28.8 

- 14.3 

16.7 

+ 24.5 

- 11.6 

4.0-4.5 

19.9 

+ 26.0 

- 14.3 

18.9 

+ 25.0 

- 13.2 

16.6 

+ 21.6 

- 11.5 

14.2 

+ 18.7 

- 9.6 

11.8 

+ 15.9 
- 7.9 

4.5-5.0 

15.2 

+ 19.0 
- 10.5 

13.9 

+ 17.3 
- 9.4 

12.1 

+ 15.3 
- 8.0 

10.3 

+ 13.0 

- 6.7 

8.3 

+ 10.9 
- 5.3 

5.0-6.0 

10.2 

+ 11.7 
- 6.7 

9.0 

+ 10.3 
- 5.8 

8.1 

+ 9.3 

- 5.0 

6.7 

+ 8.0 

- 4.1 

5.0 

+ 6.0 

- 3.0 

6.0-7.0 

5.85 

+ 6.02 

- 3.53 

5.11 

+ 5.41 

- 3.0 

4.43 

+ 4.62 

- 2.55 

3.55 

+ 3.91 

- 1.99 

2.67 

+ 2.88 

- 1.48 

7.0-8.0 

3.63 

+ 3.41 

- 2.0 

3.11 

+ 3.02 

- 1.7 

2.65 

+ 2.61 

- 1.38 

2.06 

+ 2.0 

- 1.08 

1.52 

+ 1.55 
- 0.78 

8.0-9.0 

2.21 

+ 1.98 
- 1.14 

1.89 

+ 1.68 
- 0.97 

1.57 

+ 1.42 
- 0.77 

1.25 

+ 1.2 

- 0.6 

0.78 

+ 0.74 

- 0.37 

9.0- 10.0 

1.45 

+ 1.23 
- 0.71 

1.21 

+ 1.05 
- 0.57 

1.0 

+ 0.87 

- 0.47 

0.75 

+ 0.67 

- 0.34 

0.46 

+ 0.41 

- 0.2 

10.0-15.0 

0.54 

+ 0.42 

- 0.24 

0.45 

+ 0.35 

- 0.19 

0.34 

+ 0.26 

- 0.14 

0.25 

+ 0.2 

- 0.1 

0.16 

+ 0.12 

- 0.06 

15.0-20.0 

0.12 

+ 0.08 

- 0.04 

0.1 

+ 0.07 

- 0.04 

0.07 

+ 0.05 

- 0.03 

0.04 

+ 0.03 

- 0.01 

0.02 

+ 0.02 

- 0.01 

20.0 - 30.0 

0.029 

+ 0.018 

- 0.008 

0.025 

+ 0.014 

- 0.008 

0.016 

+ 0.009 

- 0.005 

0.009 

+ 0.006 

- 0.003 

0.003 

+ 0.002 

- 0.001 


Table 3: Predictions for the differential cross-sections for meson production at LHCb at 13 TeV, 
computed using POWHEG and NNPDF3.0-|-LHCb. For each bin we indicate the central value and 
the total theoretical uncertainty. Predictions for different binnings and for other D meson species are 
available upon request. 

motivation, in Table we show ratio i?i 3 / 7 (orig) computed with the original NNPDF3.0 PDF 
set, which should be compared with the predictions obtained with the NNPDF3.0-|-LHCb set in 
Table The data is ordered in increasing rapidity bins, and within each of these in increasing 
Pt bins. For each bin, we show the central prediction, the PDF uncertainty and the the total 
theory uncertainty for i?^yy(orig), as well as the ratio between the predictions for the ratio itself 
computed with NNPDF3.0-|-LHCb and with the original NNPDF3.0, i?^y^(new)/i?^y^(orig). 

From the comparison between Tables and we see first of all that the predictions for the 
central value of are reasonably stable: differences for the central value computed between 

the original and new PDFs are typically a few percent, rather smaller than the total theory 
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p¥ (GeV) 





y 

B 






2.0 

- 2.5 

2.5 

- 3.0 

3.0 

-3.5 

3.5 

-4.0 

4.0- 

4.5 

0.0- 0.5 

0.6 

+ 0.37 

- 0.28 

0.53 

+ 0.34 

- 0.25 

0.45 

+ 0.29 

- 0.21 

0.4 

+ 0.25 

- 0.18 

0.32 

+ 0.2 

- 0.13 

0.5- 1.0 

1.67 

+ 1.0 

- 0.8 

1.53 

+ 0.94 

- 0.73 

1.36 

+ 0.84 

- 0.62 

1.14 

+ 0.71 

- 0.51 

0.91 

+ 0.56 

- 0.4 

1.0- 1.5 

2.56 

+ 1.58 

- 1.21 

2.32 

+ 1.42 
- 1.07 

2.05 

+ 1.27 
- 0.93 

1.75 

+ 1.07 
- 0.77 

1.39 

+ 0.83 

- 0.59 

1.5-2.0 

3.27 

+ 1.98 

- 1.54 

2.96 

+ 1.77 
- 1.36 

2.58 

+ 1.57 
- 1.15 

2.21 

+ 1.33 
- 0.97 

1.77 

+ 1.07 
- 0.76 

2.0- 2.5 

3.68 

+ 2.19 

- 1.68 

3.34 

+ 1.99 
- 1.49 

2.93 

+ 1.76 
- 1.29 

2.45 

+ 1.46 
- 1.04 

1.98 

+ 1.18 
- 0.83 

2.5-3.0 

3.91 

+ 2.3 

- 1.76 

3.55 

+ 2.09 

- 1.57 

3.1 

+ 1.82 
- 1.33 

2.58 

+ 1.52 
- 1.07 

2.05 

+ 1.2 
- 0.84 

3.0- 3.5 

3.99 

+ 2.25 

- 1.76 

3.58 

+ 2.06 

- 1.53 

3.16 

+ 1.82 
- 1.32 

2.64 

+ 1.51 
- 1.09 

2.04 

+ 1.16 
- 0.82 

3.5-4.0 

3.91 

+ 2.17 

- 1.68 

3.49 

+ 1.96 
- 1.46 

3.06 

+ 1.7 
- 1.25 

2.54 

+ 1.44 
- 1.01 

1.97 

+ 1.11 
- 0.76 

4.0-4.5 

3.71 

+ 2.03 

- 1.54 

3.36 

+ 1.83 
- 1.37 

2.9 

+ 1.6 
- 1.17 

2.4 

+ 1.3 
- 0.93 

1.82 

+ 0.98 

- 0.67 

4.5-5.0 

3.5 

+ 1.87 

- 1.43 

3.16 

+ 1.68 
- 1.27 

2.67 

+ 1.42 
- 1.04 

2.21 

+ 1.18 
- 0.84 

1.69 

+ 0.91 

- 0.62 

5.0- 5.5 

3.17 

+ 1.69 

- 1.25 

2.87 

+ 1.49 
- 1.12 

2.43 

+ 1.27 
- 0.91 

1.98 

+ 1.05 
- 0.72 

1.5 

+ 0.77 

- 0.53 

5.5-6.0 

2.89 

+ 1.48 

- 1.12 

2.6 

+ 1.31 
- 0.98 

2.22 

+ 1.13 
- 0.82 

1.77 

+ 0.89 

- 0.62 

1.34 

+ 0.68 

- 0.46 

6.0- 6.5 

2.63 

+ 1.29 

- 1.0 

2.33 

+ 1.15 
- 0.87 

1.97 

+ 0.97 

- 0.69 

1.58 

+ 0.76 

- 0.55 

1.15 

+ 0.57 

- 0.37 

6.5-7.0 

2.36 

+ 1.14 

- 0.87 

2.05 

+ 1.0 
- 0.73 

1.74 

+ 0.85 

- 0.61 

1.38 

+ 0.67 

- 0.45 

0.98 

+ 0.48 

- 0.31 

1 

O 

2.07 

+ 0.99 

- 0.74 

1.84 

+ 0.88 

- 0.65 

1.52 

+ 0.73 

- 0.5 

1.2 

+ 0.56 

- 0.38 

0.86 

+ 0.41 

- 0.27 

q 

00 

1 

1.82 

+ 0.87 

- 0.63 

1.61 

+ 0.76 

- 0.54 

1.33 

+ 0.62 

- 0.43 

1.04 

+ 0.47 

- 0.32 

0.72 

+ 0.34 

- 0.22 

rO 

00 

1 

q 

00 

1.62 

+ 0.74 

- 0.55 

1.41 

+ 0.65 

- 0.47 

1.15 

+ 0.51 

- 0.36 

0.88 

+ 0.39 

- 0.26 

0.62 

+ 0.28 

- 0.18 

8.5-9.0 

1.41 

+ 0.63 

- 0.46 

1.23 

+ 0.55 

- 0.39 

1.01 

+ 0.45 

- 0.31 

0.77 

+ 0.35 

- 0.23 

0.52 

+ 0.24 

- 0.14 

9.0- 9.5 

1.22 

+ 0.55 

- 0.39 

1.07 

+ 0.48 

- 0.33 

0.87 

+ 0.38 

- 0.26 

0.66 

+ 0.29 

- 0.19 

0.45 

+ 0.2 

- 0.12 

9.5-10.0 

1.11 

+ 0.48 

- 0.35 

0.93 

+ 0.4 

- 0.28 

0.77 

+ 0.33 

- 0.23 

0.58 

+ 0.26 

- 0.16 

0.39 

+ 0.16 

- 0.11 

10.0- 10.5 

0.95 

+ 0.41 

- 0.29 

0.82 

+ 0.35 

- 0.24 

0.66 

+ 0.28 

- 0.19 

0.48 

+ 0.21 

- 0.12 

0.33 

+ 0.14 

- 0.09 

10.5- 11.5 

0.79 

+ 0.33 

- 0.23 

0.67 

+ 0.28 

- 0.19 

0.53 

+ 0.22 

- 0.14 

0.4 

+ 0.16 

- 0.1 

0.26 

+ 0.11 

- 0.07 

11.5 - 12.5 

0.61 

+ 0.25 

- 0.17 

0.51 

+ 0.21 

- 0.14 

0.4 

+ 0.17 

- 0.1 

0.3 

+ 0.12 

- 0.08 

0.18 

+ 0.07 

- 0.05 

12.5 - 14.0 

0.45 

+ 0.18 

- 0.12 

0.38 

+ 0.15 

- 0.1 

0.29 

+ 0.12 

- 0.08 

0.21 

+ 0.08 

- 0.05 

0.13 

+ 0.05 

- 0.03 

14.0- 16.5 

0.28 

+ 0.11 

- 0.07 

0.23 

+ 0.09 

- 0.06 

0.18 

+ 0.07 

- 0.04 

0.12 

+ 0.05 

- 0.03 

0.07 

+ 0.03 

- 0.02 

16.5-23.5 

0.11 

+ 0.04 

- 0.03 

0.09 

+ 0.03 

- 0.02 

0.06 

+ 0.02 

- 0.01 

0.04 

+ 0.01 

- 0.01 

0.02 

+ 0.01 

- 0.0 

23.5 - 40.0 

0.019 

+ 0.006 

- 0.004 

0.014 

+ 0.005 

- 0.003 

0.01 

+ 0.003 

- 0.002 

0.005 

+ 0.002 

- 0.001 

0.002 

+ 0.001 

- 0.001 


Table 4: Same as Table now for the production of i?° mesons at 13 TeV. Predictions for different 
binnings and for other B mesons are available upon request. 

uncertainties. This nicely illustrates the compatibility of the NNPDF3.0 small-a: gluon with the 
7 TeV LHCb charm production data. The real difference comes from the reduction in PDF 
uncertainties: since scale and charm mass uncertainties are essentially the same in ii^^.j,(new) 
and i?^^y(old), the differences between the total theory errors stem from the reduction of PDF 
uncertainties in i?i 3 / 7 (new). For instance, in the lowest pT and most forward region (data 
bin 33), the relative total theory uncertainty of i?^^.j,(orig) is while for i?^^^(new) the 
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rD _ dV(D)(y,pr,13 TeV) j d^a(D)iy,pT,7 TeV) 
13/7 dy^dpi^ / dy^dp!^ 


Pt 

(GeV) 





y 

D 






2 . 0 - 

- 2.5 

2.5 

- 3.0 

3.0 

- 3.5 

3.5 

-4.0 

4.0- 

-4.5 

0.0 

- 1.0 

1.23 

-1-0.13 

-0.35 

1.25 

-1-0.13 

-0.35 

1.28 

-1-0.15 

-0.32 

1.33 

-1-0.15 

-0.32 

1.38 

+0.17 

-0.27 

1.0 

- 2.0 

1.26 

-1-0.12 

-0.34 

1.29 

-1-0.12 

-0.37 

1.32 

-1-0.13 

-0.35 

1.38 

-1-0.13 

-0.31 

1.46 

+0.14 

-0.31 

2.0 

- 3.0 

1.31 

-1-0.11 

-0.31 

1.33 

-1-0.12 

-0.29 

1.39 

-1-0.12 

-0.31 

1.45 

-1-0.13 

-0.33 

1.57 

+0.15 

-0.3 

3.0 

-4.0 

1.39 

-ro.i 

-0.26 

1.43 

-1-0.11 

-0.26 

1.5 

-1-0.11 

-0.27 

1.61 

-1-0.13 

-0.33 

1.72 

+0.16 

-0.32 

4.0 

- 5.0 

1.49 

-ro.i 

-0.22 

1.53 

-1-0.13 

-0.24 

1.61 

-1-0.12 

-0.26 

1.72 

-1-0.14 

-0.29 

2.01 

+0.16 

-0.39 

5.0 

- 6.0 

1.57 

-1-0.11 

-0.24 

1.65 

-1-0.12 

-0.25 

1.79 

-1-0.11 

-0.26 

1.99 

-1-0.11 

-0.25 

2.18 

+0.19 

-0.35 

6.0 

- 7.0 

1.67 

-1-0.13 

-0.24 

1.75 

-1-0.12 

-0.23 

1.85 

-1-0.14 

-0.28 

1.98 

-1-0.21 

-0.27 



7.0 

- 8.0 

1.78 

-ro.i 

-0.18 

1.84 

-1-0.09 

-0.18 

2.02 

-1-0.12 

-0.15 

2.17 

-1-0.16 

-0.31 




Table 5: Predictions for the ratio of double differential cross-sections for meson produc¬ 

tion between 13 and 7 TeV at LHCb, Eq. Q Results have obtained using POWHEG with the 
NNPDF3.0-|-LHCb NLO PDF set. The same binning as in the 7 TeV measurement is assumed at 
13 TeV. In each bin, we provide the central prediction and the total theoretical uncertainty, obtained 
from the sum in quadrature of scales, PDFs and charm mass variations. See Fig. for the graphical 
representation of these predictions. 


corresponding uncertainty is substantially reduced reduced down to ^ 27 % ■ Similar comparisons 
can be performed for other bins. 

We should mention that, once a measurement of R 13/7 becomes available, it should be 
possible to include this data in a global PDF fit in a similar way as we have done with the 7 
TeV charm normalised cross-sections. One expects similar improvements in the low-x gluon, 
though perhaps the increased lever arm in x of the 13 TeV data will increase the constraining 
power towards smaller values of x. As discussed in Sect. |4.2[ the main advantage of the ratio 
measurement is the cancellation of theory systematics, in particular from scale variations. 

We have also computed the value of ratio of inclusive fiducial cross-sections, as explained in 


Sect. 4.3, but this time for original NNPDF3.0 set, which turns out to be 


o f ■ \ 1 (-0+0-19(12.6%) 

7?D(orig) = 


(19) 


where we provide the total theoretical uncertainty of the prediction. This should be compared 
with the result obtained with the NNPDF3.0-|-LHCb set, Eq. Q. The reduction of the total 
theory uncertainty in Eq. ([^ as compared to Eq. (19) is a consequence of the constraints from 
the 7 TeV LHCb charm measurements. 

We now provide the differential predictions for B meson production at LHC Run H in LHCb. 
Eirst of all, in Tablej^we provide the predictions for the ratio of double differential cross-sections 
for the production of mesons between 13 TeV and 7 TeV, see Eq. ([^. These results were 
represented graphically in Fig. [TB} This is the analog table as that for charm production in 
Table As in the case of charm, we have assumed the same binning in than the 

corresponding 7 TeV measurement. The magnitude of R ^/7 increases rapidly with increasing 

and where the 7 TeV cross-sections are close to their kinematical boundaries. 
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Data Index 

i?i3/7(new)/i?i3/7(orig) 

^i3/7(orig) cv 

-Ri3/7(orig) PDF 

-Ri3/7(orig) Tot 

1 

0.9 

1.36 

±0.19 

+0.23 

-0.43 

2 

0.93 

1.35 

±0.18 

+0.22 

-0.4 

3 

0.94 

1.39 

±0.19 

+0.21 

-0.37 

4 

0.95 

1.47 

±0.18 

+0.2 

-0.32 

5 

0.97 

1.55 

±0.15 

+0.18 

-0.27 

6 

0.99 

1.59 

±0.13 

+0.17 

-0.27 

7 

1.0 

1.67 

±0.12 

+0.18 

-0.27 

8 

1.01 

1.75 

±0.09 

+0.13 

-0.2 

9 

0.88 

1.42 

±0.19 

+0.23 

-0.43 

10 

0.92 

1.41 

±0.17 

+0.21 

-0.44 

11 

0.93 

1.42 

±0.17 

+0.2 

-0.35 

12 

0.96 

1.5 

±0.16 

+0.19 

-0.31 

13 

0.96 

1.6 

±0.15 

+0.2 

-0.29 

14 

0.95 

1.73 

±0.18 

+0.22 

-0.31 

15 

0.96 

1.81 

±0.16 

+0.2 

-0.29 

16 

0.99 

1.85 

±0.14 

+0.16 

-0.22 

17 

0.88 

1.46 

±0.2 

+0.25 

-0.4 

18 

0.91 

1.44 

±0.18 

+0.22 

-0.42 

19 

0.93 

1.48 

±0.16 

+0.2 

-0.36 

20 

0.95 

1.58 

±0.14 

+0.18 

-0.32 

21 

0.95 

1.69 

±0.15 

+0.19 

-0.3 

22 

1.02 

1.75 

±0.16 

+0.19 

-0.3 

23 

1.0 

1.85 

±0.16 

+0.21 

-0.32 

24 

1.02 

1.97 

±0.17 

+0.19 

-0.21 

25 

0.88 

1.51 

±0.25 

+0.29 

-0.43 

26 

0.91 

1.51 

±0.19 

+0.23 

-0.38 

27 

0.92 

1.57 

±0.16 

+0.21 

-0.39 

28 

0.96 

1.68 

±0.15 

+0.2 

-0.37 

29 

0.97 

1.78 

±0.15 

+0.2 

-0.33 

30 

1.01 

1.98 

±0.16 

+0.18 

-0.29 

31 

0.95 

2.09 

±0.17 

+0.27 

-0.32 

32 

1.04 

2.09 

±0.16 

+0.21 

-0.33 

33 

0.86 

1.61 

±0.24 

+0.3 

-0.38 

34 

0.9 

1.61 

±0.21 

+0.25 

-0.39 

35 

0.92 

1.7 

±0.18 

+0.23 

-0.36 

36 

0.93 

1.85 

±0.17 

+0.23 

-0.37 

37 

0.98 

2.04 

±0.16 

+0.22 

-0.42 

38 

0.99 

2.21 

±0.2 

+0.27 

-0.4 


Table 6: The ratio R^ 3 /^ computed with the original NNPDF3.0 PDF set, in order to compare with 
the predictions obtained with the NNPDF3.0+LHCb set in Table The data is ordered in increasing 
rapidity bins, an within each of these, in increasing pT bins. For each bin, we show the central prediction, 
the PDF uncertainty and the the total theory uncertainty for iii 3 / 7 (orig), as well as the ratio between 
the new and orig predictions for the ratio itself, i?i 3 / 7 (new)/i?i 3 / 7 (orig). 
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oS _ d'^ <7 (B)[y,pT,13 TeV) / d? tT(B){y ,PT ,7 TeV) 
1317 _ dyB dpB/dyB dpB 


p¥ (GeV) 



2.0 

- 2.5 

2.5 
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3.0- 

-3.5 

3.5 

-4.0 

4.0- 

4.5 
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1.6 

+ 0.08 
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+ 0.08 

- 0.14 
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1.98 

+ 0.13 
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3.0-3.5 
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+ 0.07 

- 0.13 
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- 0.14 
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- 0.14 
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- 0.13 

1.65 

+ 0.07 

- 0.14 
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- 0.15 
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+ 0.1 
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- 0.19 

7.0-7.5 

1.74 

+ 0.07 

- 0.14 
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1.95 
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2.23 

+ 0.08 

- 0.17 

2.54 
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- 0.21 

7.5-8.0 

1.76 
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- 0.13 

1.85 
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1.97 
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- 0.14 

2.24 
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- 0.17 

2.53 

+ 0.19 

- 0.2 

8.0-8.5 

1.81 
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1.93 
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2.03 
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- 0.14 
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- 0.21 
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